TWO-DIMENSIONAL SUPERINTEGRABLE METRICS WITH 
ONE LINEAR AND ONE CUBIC INTEGRAL. 



VLADIMIR S. MATVEEV AND VSEVOLOD V. SHEVCHISHIN 

Abstract. We describe all local Riemannian metrics on surfaces whose geodesic flows are 
\ superintegrable with one integral linear in momenta and one integral cubic in momenta. 

We also show that some of these metrics can be extended to S 2 . This gives us new 
examples of Hamiltonian systems on the sphere with integrals of degree three in momenta, 
■ and the first examples of superintegrable metrics of nonconstant curvature on a closed 

\ surface. 
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1. Introduction 



43 ■ 1.1. Definitions and statement of the problem. Let M 2 be a surface (i.e., 2- dimensional 
real manifold) equipped with a Riemannian metric g = (gij). The geodesic flow of the metric 
g is the Hamiltonian system on the cotangent bundle T*M 2 with the Hamiltonian H : = 
\g l ipiPj, where (x,y) = (xi,X2) is a local coordinate system on M 2 , and (p x ,p y ) = (pi,P2) 
are the correspondent momenta, i.e., the dual coordinates on T*M 2 . 

We say that a function F : T*M 2 — > M. is an integral of the geodesic flow of g, if {F, H} = 0, 



where { , } is the canonical Poisson bracket on T*M 2 . We say that the integral is polynomial 
in momenta of degree d, if in every local coordinate system (x,y,p x ,p y ) it has the form 

F(x,y,p x ,p y ) = Y^=o a *( x > v)P%~ i P i y i 
For example, the Hamiltonian H itself is an integral quadratic in momenta. Integrals 
polynomial in momenta of degree 1 (3, resp.) will be called linear (resp. cubic) integrals. 
The first main result of the present paper is a complete solution of the following problem: 

X " 

5^ ■ Problem A. Describe locally all two-dimensional Riemannian metrics admitting one inte- 
gral L linear in momenta and one integral F cubic in momenta such that L, F, and H are 
functionally independent. 

Recall that functions L, F, H are functionally independent if there exists a point on T*M 
such that at this point the differentials dL, dF, dH are linearly independent. For integrals 
polynomial in momenta, linear independence of the differentials of the integrals at one point 
implies linear independence of the differentials of the integrals at every point of a certain 
everywhere dense open subset (assuming the manifold is connected). 

Recall that two-dimensional metrics whose geodesic flows admit three functionally inde- 
pendent integrals of a certain special form (in most cases the integrals are assumed to be 
polynomial in momenta of certain fixed degrees) are called superintegrable. Superintegrable 
metrics (and Hamiltonian systems in general) is nowadays a hot topic in mathematical 
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physics and differential geometry, due to various applications and interesting mathematical 
structure lying behind. We suggest |K-K-M-W1 IK-K-M} ITe-Wi} |Da-Yp| for a discussion of 



superintegrable systems from the viewpoint of mathematical physics, and \Kt\ IBr-Ma-Ma] 
from the viewpoint of differential geometry. 

If the metric is superintegrable with two linear integrals, it has constant curvature. The 
metrics that are superintegrable with two quadratic integrals (in addition to the energy 
integral), or one linear and one quadratic, were described (locally, in a neighbourhood of 
almost every point) in the classical work of Koenigs |Koe] . 

The next case should be "linear integral + cubic integral" , but the only result we found in 
this direction is due to Ranada |Raj . Gravel |Grj . and Marquette and Winternitz |Mar-Wi] 
and concerns the Hamiltonian systems such that the Hamiltonian H is the sum of the stan- 
dard kinetic energy K = \{p 2 x + p 2 ) and a potential energy V(x, y) ^ const. They assumed 
the existence of (functionally independent) linear and cubic integrals and proved that for 
such systems the cubic integral is actually the product of the linear integral and of an integral 
quadratic in momenta, i.e., such systems can be obtained via the Maupertuis' transformation 
from the superintegrable systems constructed by Koenigs [Koej . In particular, all known ex- 
amples of metrics satisfying assumptions in Problem A above were in a certain sense trivial: 
the metric has a constant curvature, or the metric is superintegrable with one quadratic and 
one linear integral, and every cubic integral is a product of the integral linear in momenta 
and an integral quadratic in momenta. 

1.2. Main result: local normal forms of metrics admitting one linear and one 
cubic integral. 

Theorem 1.1. Let g be a Riemannian metric on the 2- dimensional connected manifold M 2 . 

Suppose the geodesic flow of g admits a linear integral L and a cubic integral F such that 
L, F and the Hamiltonian H are functionally independent. Then, locally near every point 
p such that L\ T * M 2 ^ there exist coordinates (x,y) and a real function h(x) such that the 
metric g has the form g = -^-(dx 2 + dy 2 ) where h x is the x-derivative, and the function 
satisfies one of the following Principal (ordinary differential) equations: 

(i) h x -(A -h 2 x + fi 2 -A -h{x) 2 -A v h{x) + A 2 ) - (A 3 - 5 -^ + JU-cosQjl-x)) =0 

(it) h x -(A Q -h 2 x - ^-A -h{xf - A r h{x) + A 2 ) - (A 3 - s -^f^ + A 4 -cosh(/i-x)) =0 (1.2) 

(Hi) h x -(A -h 2 x - A v h(x) + A 2 ) - (A 3 -x + A 4 ) = 

with \i > in the first two cases. 

In all three cases the metric g = -^(dx 2 + dy 2 ) is superintegrable with one linear integral 
L = p y and one cubic integral. In the case (i) a cubic integral F = F(x,y;p x ,p y ) can be 
given by 

F = (C + -e^ + C_ - e-«0 • (a (x) -pi + a 2 (x) -p x p 2 y ) + (C + -e™ - C_ -e^) • (a x (x) -p 2 x p y + a 3 (x) -pf) 
where C + ,C_ are arbitrary constants and ai(x) are functions given by 
a {x) = A -h 3 x 



ai {x) = {-^A -h{x) + ^)-h 2 x 



a 2 (x) = \ i3A -h 2 x + fi 2 -A -h(x) 2 - A v h(x) + A 2 )-h x 
a 3 (x) = ±-{3A -h 2 x + fi 2 -A -h(x) 2 - A v h(x) + A 2 )-h x 



1.3) 
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In the case (ii) a cubic integral can be given by 

F = C e -cos(fi-y + (j) )-(a (x)-pl + a 2 (x)-p x p 2 y ) + C e -s\n(fi-y + <j) )-{ai(x)-plp y + a 3 (x)-p 3 y ) 
where C e ,<f)o are constants ("amplitude and phase") and a,i(x) are functions given by 
a (x) = A -h 3 x 

ai(x) = {(i-Ao-h(x) + ^)-h 2 x 

a 2 (x) = \ -{3Ao-h 2 x - n 2 -A -h(x) 2 - A v h(x) + A 2 )-h x 
a ^ x ) = Jili 3A o-h 2 x - fi 2 -A -h(x) 2 - A r h{x) + A 2 )-h xx . 
In the case (Hi) a cubic integral can be given by 

F = C v (a (x)-pl + a 2 (x)-p x p 2 y + §• (A v h 2 x -p 2 xPy + (A v h 2 x + A 3 )-p 3 y ) 

+ G T (y-a (x)-pl + ai(x)-j^p w + y-a 2 (x)-p x p 2 y + a 3 (x)-p 3 y (1.5) 

+ y i\M-h 2 x -p 2 x p y + {Ayh 2 + A 3 )-p 3 y ) 

where C±, C 2 are constants and ai(x) are functions given by 

a (x) = A -h 3 x 
ai(z) = -A -h 2 x -h(x) 
a 2 (x) = l-{3-h 2 x -A - A v h(x) + A 2 )-h x 
aa(x) = -\-{4:A -h 2 x -h{x) + A 3 -x 2 + 2A r x). 

Moreover, in the case when the metric g has non-constant curvature every cubic integral is 
a linear combination F+Cl 3 -L 3 + Clh-L-H where F is given by the above formula (according 
to the case (i)-(iii)) and Cl 3 ,Clh o,re constants. In particular, in the non-constant curvature 
case the space of cubic integrals of our metrics has dimension 4. 

Remark 1.1. Uniqueness of the equation. We show in Theorem 14.11 that in the case 
when the curvature of our metric g = -^(dx 2 + dy 2 ) is non-constant the equation ( II. 2ft on 
the function h(x) is unique up to a constant factor. On the other hand, in Theorem 13.11 we 
describe possible equations of the form (11.21) for whose the metric g = -^(dx 2 + dy 2 ) has 
constant curvature. 

Thus, Theorems 11.11 13.11 and 14.11 give a complete answer to the Problem A above. 

Remark 1.2. Known special case: Darboux-superintegrable metrics. We call a 
metric g on M 2 Darboux-superintegrable, if it has non-constant curvature and the geodesic 
flow of the metric admits at least four linear independent integrals quadratic in momenta. 
For example, such is the metric (x 2 + y 2 + l)-{dx 2 + dy 2 ) on M 2 (see for example |Malt §4]). 

Darboux-superintegrable metrics are well understood. Locally, they were described already 
by Koenigs |Koe] . In particular, Koenigs has shown that every Darboux-superintegrable 
metric admits a linear integral. Then, it also admits cubic integrals, namely the products 
of the linear integral and the quadratic integrals. Therefore, our Theorem 11.11 applies. In 
particular, in the appropriate local coordinates (x, y) the metric has the form -^{dx 2 + dy 2 ) 
such that h(x) is a solution of the Principle equation (jl.2p . (i)-(iii). 

The formulas above show that if the coefficient Aq vanishes, then a generic cubic integral 
F + Cl 3 -L 3 + Clh-L-H is the product of the linear integral and a function quadratic in 
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momenta which must automatically be an integral. Then, the metrics corresponding to 
A = are Darboux-superintegrable or are of constant curvature. The uniqueness of the 
Principle equation (see the previous remark) shows that the converse assertion is also true. 
Thus we obtain the following characterisation: Under the hypotheses of Theorem li.il the 
metric g = -^(dx 2 + dy 2 ) of non-constant curvature is Darboux-superintegrable if and only 
if the parameter A vanishes. 

Remark 1.3. With the help of a computer algebra software, for example with Maple®, it 
is easy to check that the functions F from Theorem 11.11 are indeed integrals: the condition 
{H, F} = is equivalent to 5 ODEs of at most 3rd order on the function h; these ODEs are 
identically fulfilled for the function h satisfying the corresponding equation (jl.2p . since they 
are algebraic corollaries of the corresponding equation and its first two derivatives. We will 
of course explain how we constructed the integrals (in Section [2]) since we also need to show 
that we constructed all such metrics. Moreover, the idea of the construction will be used in 
the proof of other statements of Theorem 11.11 in particular in the proof that the dimension 
of the space of the integrals is 4. Moreover, we believe that the idea of our construction 
could also be used for constructing higher order superintegrable cases, see Problem [1] in the 
conclusion. 

1.3. Second main result: Examples of metrics on the 2-sphere admitting linear 
and cubic integrals. The problem of finding and describing global integrable Hamiltonian 
systems, i.e., those whose configuration space is a compact manifold, is one of the central 
topics in the classical mechanics. The version of this problem in our context is as follows: 

Problem B. Understand what Riemannian metrics on the 2-sphere S 2 admit one integral L 
linear in momenta and one integral F cubic in momenta such that L, F, and the Hamiltonian 
H are functionally independent. 

Note that other oriented closed surfaces can not admit superintegrable metrics. Indeed, 
if the metric is superintegrable, all geodesies are closed which is possible on the sphere and 
on RP 2 only. 

It is known that the metric of constant curvature of the sphere do admit (linearly inde- 
pendent) linear and cubic integrals. So the nontrivial part of the Problem B is whether there 
are other metrics on the 2-sphere admitting an integral L linear in momenta and an integral 
F cubic in momenta such that L, F, and H are functionally independent. 

For the integrals of lower degrees, the answer in negative. Indeed, the existence of two 
functionally independent linear integrals implies, even locally, that the metric has constant 
curvature. By Kiyohara |Kij . the existence of three functionally independent quadratic in- 
tegral (energy integral + two additional integrals) on the 2-sphere implies that the metric 
is of constant curvature. From this result, it also follows that the existence of (function- 
ally independent) linear and two quadratic integrals implies that the metric of the 2-sphere 
has constant curvature. Because of these results (and absence of examples of superinte- 
grable systems with higher degree integrals), it was generally believed that no polynomially 
superintegrable metric exists on a closed surface of nonconstant curvature. 

In the present paper, we construct the first examples of smooth (even analytic) metrics 
of nonconstant curvature on the 2-sphere whose geodesic flows admit integrals L linear in 
momenta and integrals F cubic in momenta such that L, F, and H are functionally indepen- 
dent. The construction is in Section [61 We show that for certain values of parameters the 
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metrics we constructed in Theorem 11.11 and the integrals of these metrics can be smoothly 
extended to the sphere. More precisely (in the notation of Theorem ll.ip . if the function 
h(x) fulfills the equation ( II. 2p (ii) and the condition h'(xo) > at some point xq whereas the 
real parameters ji > 0,A , . . . , At satisfy inequalities A > 0, [i-A^ > \A 3 \ then the metric 
g = j2 (dx 2 + dy 2 ) smoothly extends to the sphere S 2 together with the linear integral L = p y 
and the cubic integral F given by (II .4p . 

We conject that these are all examples of metrics on the sphere superintegrable with one 
linear and one cubic integral. 

Our examples are also interesting from other points of view. Indeed, every metric from 
these examples admit and integral cubic in momenta that not the product of a linear and 
a quadratic integral. The problem of constructing such metrics is very classical and, in 
a certain extend, was stated by Jacobi, Darboux, Cauchy, Whittaker, see also |Bo-Ko-Fo"l 
IBo-Foj . There are only very few examples of such metrics on closed surfaces: constant 
curvature metrics, metrics constructed via Maupertuis' transformation from the Goryachev- 
Chaplygin case of rigid body motion and their generalizations due to Goryachev |Go] , metrics 
constructed by Selivanova j5e] , by Kiyohara |Ki2j . and by Dullin at al |Du-Maj (see also 
|Du-Ma-Toj and |Ma-Shj ). Note that the analo gous question for the quadratic integrals is 
completely solved, see [KoH EH IBo-Ma-Fol IMaj . 

Moreover, all geodesies of the metrics we constructed are closed (since it is always the case 
for superintegrable metrics), so the examples are also examples of the so-called Zoll surfaces. 

1.4. Additional result: Special case of Kruglikov's "big gap" conjecture. In |Kr| 

§12], Kruglikov has shown that the dimension of the space of cubic integrals (of the geodesic 
flow of a 2D-metric) is at most 10; the dimension 10 is achieved only by the metrics of 
constant curvature. He also has shown that the second largest dimension is at most 7 (see 
|Kr| Theorem 8]), and conjectured that the gap between the largest and the second largest 
possibilities for the dimension of the space of cubic integrals is even bigger: he writes that 
it seems that the next realized dimension after 10 is 4- 

We will prove this conjecture (see Theorem 15. ip under the additional assumption that the 
metric admits a Killing vector field. Note that this assumption does not look too artificial, 
since it is expected that metrics with many polynomial integrals admit Killing vector fields. 
For example, by the classical result of Koenigs mentioned above, metrics admitting four (= 
the second largest dimension) linearly independent integrals that are quadratic in momenta 
admit Killing vector fields. 

More precisely, we will prove that, if a 2D metric is superintegrable with one linear and one 
cubic integral (L and F) and has non-constant curvature, then, locally, the space of cubic 
integrals is precisely 4-dimensional. In particular, in addition to the integrals L 3 , F, L ■ H we 
always construct one more cubic integral F 2 that is linearly independent of L 3 , F, and L ■ H . 

2. Principle equation and overview of the proof of Theorem 11.11 

2.1. How we found the metrics: scheme of the proof of Theorem 11.11 It is well- 
known (see for example |Da[ §592], or |Bo-Ma-Fo] ) that every pair (g,L), where g is a Rie- 
mannian metric, and L is an integral linear in momenta, is given in appropriate coordinates 
in a neighbourhood of every point such that L ^ by the formulas 

g = \(x)(dx 2 + dy 2 ) and L = p y . (2.1) 
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The natural "naive" method to solve the Problem A would be to write the condition 
{H, F} = 0, where H = |^ and F := a (x, y)p 3 x + a x (x, y)p 2 x p y + a 2 (x, y)p x p 2 y + a 3 (x, y)p 3 y , 
as the systems of PDE on the unknown function A of one variable and unknown functions a, 
of two variables, and to try to solve it. Unfortunately, by this method we obtain a system of 
5 nonlinear PDE on 5 unknown functions A, do, a±, a 2l 03, which is completely intractable^. 

In order to solve the problem, we used a trick that allowed us to reduce the problem to 
solving systems of ODE (instead of PDE). A similar trick was recently used in |Ma2] . 

The main observation is the following: the Poisson bracket of the linear integral L and of 
a cubic integral F is 

• an integral (because of the Jacobi identity) , and 

• is cubic in momenta (because each term in the sum {L, F} = d x Fd Px L + d y Fd Py L — 
d x Ld Px F — d y Ld Py F is cubic in momenta). 

Thus, the mapping £ : F H- {L, F} is a linear homomorphism. By |Krj . the space of cubic 
integrals is finite- (at most, 10-) dimensional. Let us now consider the eigenvalues of the 
mapping C. Clearly, is an eigenvalue of C, whose eigenvectors are A 3 ■ L 3 + A\ ■ L ■ H, 
where A\, A3 £ R. The following two cases are possible: 

Case 1: The mapping C has an eigenvalue \l 7^ 0. Then, there exists a cubic integral 
F such that {L, F} = fi ■ F . We allow fi to be a complex number, and F to be complex 
valued function, i.e., F = F\ + \F 2 for real- valued cubic integrals F\ and F 2 . 
In the coordinates such that (g, L) are given by ( 12. ip . we have 

{L, F} = d y F = d y a (x, y)-p 3 x + d y ai(x, y) ■ p\p y + d y a 2 (x, y) ■ p x p\ + d y a 3 (x, y) ■ p 3 y , 

so that the condition {L, F} = fi ■ F is equivalent to the system d y ai(x,y) = fi ■ ai(x,y), 
i = 0, ...,3. Then, a,i(x,y) = exp(/xy) • aj(x) for certain (complex valued in the general case) 
functions ai{x) of one variable x. Then, all unknown functions in the equation {H, F} = 
are functions of the variable x only, i.e., the condition {H, F} = is a system of ODE 
(depending on the parameter fi). Finally, the condition {H, F} = is equivalent to 5 ODE 
on 5 unknown functions of one variable x: four unknown functions a,i(x) and A(x). Working 
with this system of ODE, we partially integrate it and reduce it to one ODE of the first 
order (essentially, the first equation of ( II. 2p for fi £ K. and the second equation of (II. 2p for 
fi £ i ■ R. Further in §4.11 we shall show that the assumption that the function h x is real 
implies that fi is real or pure imaginary). 

Remark 2.1. In this way, we obtain X(x) which is a priori a complex- valued function; for 
our problem, only real-valued A's are of interest. We shall see in §4.11 that A is real if and 
only if fi is real or purely imaginary. 

Case 2: The mapping C has only one eigenvalue, namely zero. Since in our setting 
the space of cubic integrals is at least three-dimensional, there exists an integral F linear 
independent of L 3 and L-H such that {L, F} = 4f-L 3 + Ai-L-H for certain constants A\, A3. 

1 This "naive" approach to this problem was tried without success by many experts in superintegrable 
systems (private communications by Marquette, Rahada, Winternitz). 
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In the coordinates such that (g, L) is given by flUD , the condition {L, F} = 4* ■ L 3 + A l ■ L ■ H 
reads 

{p y ,a (x,y)p 3 x + a 1 (x J y)p 2 x p y + a 2 (x,y)p x p 2 y + a 3 (x,y)p y } 
= d y a (x, y)-pl + d y a x (x, y) ■ p 2 x p y + d y a 2 (x, y) ■ p x p\ + d y a 3 (x, y) ■ p\ 



and is equivalent to the system d y ao(x, y) = 0, d y ai(x, y) = 2^y ; d y a 2 (x, y) = 0, d y a 3 (x, y) 



_l_ ai xhen 

^ 2A(x) • ^ ' 



a (z)-^ + a x {x)-p xVy + a 2 (x)-p x p 2 + a 3 (z)-pj + | ■ ^-pj + A r p y - P -^) . (2.2) 

We again see that all unknown functions in the equation F} = are functions of the 
variable x only, i.e., the condition {H, F} = is a system of 5 ODE (depending on the 
parameters Ai,A 3 ,yo) on 5 unknown functions of one variable x: at(x) and A(x). Working 
with this system of ODE, we partially integrate it and reduce it to one ODE of the first order 
(Equation (iii) of ( 11.2j) ). which is in a certain sense a degenerate case of the corresponding 
ODE we obtained in Case 1. 

2.2. Case 1 (/i ^ 0) in the proof of Theorem 11.11 For convenience in further computa- 
tion, we write the metric g in the form 

9 = (2.3) 

for some function h = h(x), where h x = Then H = ^\p\ + p 2 y ) and the linear integral 

is L := p y . 

We assume (see §2.ip that there exists a complex-valued cubic integral of the form 

F = exp(n-y)-a (x)-p 3 x + exp(/i-y)-ai(x)-p^ + exp(fi-y)-a 2 (x)-p x p 2 y + exp(^-j/)-a 3 (x)-pJ, 

where &i £1X6 smooth complex-valued functions of one real variable x. 
Then, the condition {F, H} = reads 

{F,H} = h x -exp(fj,-y)-(h x -a (x) x - 3-a Q (x)-h xx )-p x 

+ h x -exp(p-y)-(-2 ■ ax{x)-h xx + h x -(A-a (x) + h x -ai(x) x )-pl-p y 

+ h x -exp(iJ l -y)-(h x -iJ l -a 1 (x) - 3-a (x)-h xx - a 2 (x)-h xx + h x -a 2 [x) x )-p 2 x p 2 y (2.4) 

+ h x -exp(p-y)-(h x -fi-a 2 (x) + h x -a 3 (x) x - 2 • a x {x)-h xx )-p x -p z y 

+ lh x -exp(fj,-y)-(-a 2 {x)-h xx + h x -fi-a 3 (x))-p y , 

where subscripts a n (x) x ,h x mean derivation in x, and h xx is the second derivative. Since 
the monomials p x ~ l p y form a basis of homogeneous polynomials of degree 4, every line in 
( 12. 4p should vanish. This gives us a system of 5 ODEs on 5 functions h(x), a Q (x), . . . , a 3 (x): 
each line of (12.4p corresponds to one ODE. Subsequently solving the first three of them and 
resolving a 3 (x) from the last one we obtain 



a {x) = A 



A n - hi 



X 



ai(z) = (-v-A -h{x) + gfi-hl 

a 2 (x) = \-{-A v h{x) + fi 2 -A -h(x) 2 + 3-h 2 x -A + A 2 )-h x 
a 3 (x) = ±-(3-h 2 x -A - A v h(x) + p 2 -A -h(x) 2 + A 2 )-h xx , 
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with some constants A , A±, A 2 . Substituting in the remaining equation (. . .)p x -Py, we obtain 
the following non-linear ODE of order 3 on h(x): 

(3-Ao-hl + fi 2 -A -h(x) 2 - A r h{x) + A 2 )-h xxx + 
+ 6-h 2 xx -h x -A + (6-fi 2 -A -h(x) - 3-AJ-h^h^ + (2.6) 
+ 3-fi 2 -A -h 3 x + (^-A -h(x) 2 - fi 2 -A v h(x) + fjt 2 -A 2 yh x = 0. 
By direct calculations we see that the equation (12. 6p can be written in the form 

{£ + A* 2 ) {h x -(hl-A + fi 2 -A -h(x) 2 ~ A v h{x) + A 2 )) = 0. (2.7) 
Therefore this equation is equivalent to the equation 

K- (h 2 x -A + ^-AoMx) 2 - A v h(x) + A 2 ) - (A r s -^ + A 4 -cos(/ix)) = (2.8) 

in the sense that h(x) satisfies the equation (12. 6p if and only if it satisfies (12.81) with the same 
constant Aq,Ai,A 2 G C and some constants A^^A^ G C. Later, in §2.5 (see Theorem 14. 2p 
we shall show that only real Aq, A\, A 2 , A3, A4 are interesting for our purposes. 

2.3. Case 2: fi = 0. We proceed as we explained in §2.11 we write the metric in the form 

9 — jp-{.dx 2 + dy 2 ), so that now H = ^f-(p 2 +P 2 ), an d then substitute (12. 2 p in the condition 
{F, H} = 0. We obtain 

{F,H} = h x ' 

+ h x 



(h x -a (x) x - 3-a (x)-h xx )-p x 

3 



{h x -ax(x) x - 2-a 1 (x)-h xx )-p x -p y 

+ h x -{\-A r hl - 3-a (x)-h xx - a 2 (x)-h xx + h x -a 2 {x) x )-p 2 x p 2 y (2.9) 

+ h x -(h x -a 3 (x) x - 2-a 1 (x)-h xx )-p x -pl 

+ h x -(A r h 3 x + A 3 -h x - 2-a 2 {x)-h xx ))-p^ 

with the same constants A±,A 3 as in (12. 2p . This time we can subsequently resolve all 
functions ao(x), . . . , a 3 (x) from the equations and obtain 

a Q (x) = A -hl 

ai(x) = \-A r h 2 

(2 10) 

a 2 {x) = \-{3-h 2 x -A Q - A r h{x) + A 2 )-h x 
a 3 (x) = \-h 2 x -A x + A 3 , 

with some constants Aq, Ai, A 2 , A 3 . The notation in the formula above, especially A\,A 3 , 
is chosen for convenience in future formulas. Then the bracket yields 

{F, H} = -h 2 x -(3-A -h 2 x -h xx - A v h 2 x - Ayh xx -h{x) + A T h xx - A 3 )-p*. 

This means that the equation on h(x) is 

3-A -h 2 x -h xx - Ayh 2 x - A r h xx -h(x) + h xx -A 2 = A 3 . (2.11) 

The left hand side of this expression is the x-derivative of the expression 

h x -{A -hl-A r h{x)+A 2 ). (2.12) 

Therefore the equation (12. lip is equivalent to 

h x -{A -h 2 x - A r h{x) + A 2 ) - (A r x + A A ) = 0. (2.13) 
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Remark 2.2. Obviously, we obtain this equation from both the equations ( 11. 2p (i) and (ii) 
taking the limit /x — > 0. Moreover, the solution of the Cauchy initial value problem for the 
equations (12.81) . (12. 131) depends analytically on all parameters: the variable x, parameters 
/x, A , ...,A 4 , the initial point x , and the initial value h(x ). Therefore we can consider 
real solutions h(x) of the equation ( jl.2p as "real forms" of a single holomorphic multi- valued 
function h(x; \i; A , . . . , A 4 ; x , h ) depending holomorphically on the involved parameters. 
Notice also that the equation (11.21) (ii) is obtained from (ll.2p (i) by replacing \i by i -/x, and 
similarly for the corresponding cubic integrals. 

Remark 2.3. The above argumentation shows the existence of one non-trivial cubic integral 
F in the case /x = (i.e., for the metric (12 .3p with h satisfying (jl.2f iii)). Namely, such F 
can be obtained substituting the formulas (I2.10p in (12. 2p . The solution F obtained in this 
way has the form F = A 3 -L 3 + AyL-H + i*\ with a fixed cubic integral F\ which is linear 
in y. On the other hand, for \x ^ in the both cases (i) and (ii) we obtain two non- 
trivial cubic integrals linearly independent of L 3 and L-H, namely, by replacing /x by —\x in 
formulas (I2.5p . It appears that also in the case /x = there exists another cubic integral F 2 
that is (inhomogeneous) quadratic in y. The latter property is equivalent to the condition 
£ 3 (-^2) = 0. We show the existence of such F 2 in the proof of Theorem 15.11 This additional 
integral F 2 is already included in the formulas in Theorem 11.11 

This fact is the reason for the difference in formulas (II .6p and ( I2.10p . Namely, the sub- 
stitution of (11. 6p in ( 11. 5p yields the linear combination C\F\ + F 2 C 2 of two cubic integrals 
Fx, F 2 which are linear independent of L 3 and L-H. On the other hand, the substitution of 
(12717)]) in ( 12721) yields the linear combination F = A 3 -L 3 + Ax- L-H + Fx with the same Fx, 
which gives only one cubic integral linear independent of L 3 and L-H 

2.4. Remaining steps of the proof. As we have shown, if a surface metric g admits a 
linear and a non-trivial cubic integral F, then in appropriate coordinates it has the form 
h~ 2 {dx 2 + dy 2 ) for some function h(x) satisfying one of the forms (II. 2p of the Principle 
equation, and that the cubic integral F can be constructed using the formula (12. 5p or resp. 
( 12.1 Op . The remaining steps of the proof are the following: 

• We analyse in which cases the constructed metric g = h~ 2 (dx 2 + dy 2 ) belongs to already 
known types: Metrics of constant curvature and Darboux-superintegrable metrics. This 
is done in Section [3j We show that our metrics are indeed new examples for most values 
of the parameters (the values of the parameters corresponding to previously known cases 
are solutions of certain algebraic equations). 

• In Section H] we prove that in the case of non-constant curvature the function h(x) 
satisfies a unique up to constant factor equation of type (ll.2p . This result is used to 
prove the fact that if the solution h(x) of the Principle equation with complex pa- 
rameters fi; Aq, . . . , A/x is real-valued, then the parameter /x (which could be apriori 
arbitrary complex number) must be real, purely imaginary, or zero, whereas the pa- 
rameters A Q , . . . ,Ax becomes real after the application by the appropriate constant. 
This explains why we have only 3 types (i)-(iii) of the Principle equation ( II. 2p . 

• In Section \5\ we prove that in the case of non-constant curvature every type (i)-(iii) 
of Theorem 11.11 the space of cubic integrals has dimension 4. This means that un- 
der hypotheses of the main theorem there are exactly 2 non-trivial independent cubic 
integrals, in addition to L 3 and LH. This fact in a special case of Kruglikov's "big 
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gap" conjecture (see |Kr] ) about possible dimensions of the spaces of cubic integrals of 
surface metrics. 

3. Special solutions. 

In this section we consider two special cases of the Principle equation corresponding to 
Darboux-superintegrable metrics and constant curvature metrics. 

3.1. The case A = corresponds to Darboux-superintegrable metrics. Recall that 
a two-dimensional metric g is Darboux-superintegrable, if the space of its quadratic integrals 
is at least 4-dimensional and the curvature is non-constant. We shall use the following 
statement which follows from |Krj (or even from [Koej ): if a metric g (with the Hamiltonian 
H) of non- constant curvature admits a linear integral L and a quadratic integral Q such that 
L, Q and H are functionally independent, then g is Darboux-superintegrable. 

This statement implies that for every real solution h(x) of one of the equations (11.21) with 
A = the metric g = h~ 2 {dx 2 + dy 2 ) is Darboux-superintegrable. 

Indeed, Aq = if and only if the integral F from Theorem 11.11 has zero coefficient at 
p\. Since the linear integral L in Theorem 11.11 is p y , the function Q := F/p y is an integral 
quadratic in momenta. If L, H and F are functionally independent, then the functions 
L, H, Q are also functionally independent and the metric is Darboux-superintegrable by the 
result of |Koel IKrj recalled above. 

For further use let us note that every Darboux-superintegrable metric always has the form 
g = h~ 2 {dx 2 + dy 2 ) for some function h(x) satisfying one of the Principle equations ( 11. 2ft with 
A = 0. Indeed, for given metric g admitting a non- vanishing linear integral L there exists 
a isothermic coordinate system (x,y), unique up to translations, in which L = p y . In these 
coordinates g has the form g = h~ 2 (dx 2 + dy 2 ) with some function h(x). Further, if Q is a 
quadratic integral, then F := Q-L is a cubic integral for g. In this situation we have shown 
that h(x) must satisfy one of the equations (II. 2p with certain parameters A Q , . . . , A4 such 
that F = Y^H=o a i( x ^y)PxPy~ z with oq = A Q h x . The condition F = Q-L means the vanishing 
of ao(x,y) which is equivalent to A = 0. 

3.2. Parameters in Theorem 11.11 corresponding to metrics of constant curvature. 

The goal of this subsection is to understand for what values of the parameters A , v4 4 and 
the initial value h(xo) the metric from Theorem 11.11 belong to the previously known classes, 
that is to the Darboux-integrable metrics and to the metrics of constant curvature. In §3.11 
we have shown that Darboux-superintegrable metrics are characterized by the condition 
Aq = 0. Thus in order to understand whether the metrics we constructed are new we need 
to understand which metrics with Aq 7^ have constant curvature. The answer is given in 
Theorem 13.11 In particular, Corollary 13.11 shows that most metrics we constructed are new. 

Let g be a metric on M 2 of the constant Gauss curvature R and v a Killing vector 
field corresponding to the linear integral L v . Then according to the sign of R the Lie 
algebra of Killing vector fields on M 2 is either so(3) (case R>0), or s[(2,M) (case R<0), or 
the affine algebra aff(IR 2 ) of isometries of M 2 isomorphic to a semi-direct sum so (2) ix M 2 
(remaining case R = 0). The classification of elements of these three Lie algebras gives 6 
types of Killing vector fields: rotations of S 2 (R>0), rotations, hyperbolic translations, and 
loxodromies of the hyperbolic plane (R<0), and rotations and translations of M 2 (R=0). Fix 
a coordinate system (x, y) in which the metric has the form g = h~ 2 (dx 2 + dy 2 ) , the curvature 
is R = h xxx -h x — h 2 x , and the Killing vector field has the form v — Since in each of 



TWO-DIMENSIONAL METRICS WITH LINEAR AND CUBIC INTEGRALS 



11 



these cases the metric has 3 Killing vector field, there exists a cubic integral independent of 
L 3 and L v -H. Consequently, h(x) must satisfy one of the Principal equations. The explicit 
situation is as follows: 

Theorem 3.1. Assume that a metric g = h~ 2 (dx 2 + dy 2 ) has constant Gauss curvature 
R. Then h(x) satisfies one of the equations (11.21) . Moreover, in this case under additional 
assumption A ^ one of the following possibilities holds: 

(1) h(x) = a-sinh(/i-(x — b)) + c with some constants fi > 0, a > 0, c, b satisfying R = a 2 /j, 4 . 
In this case the Killing vector field is locally a rotation of the 2-sphere of radius 
r = a -1 yU -2 and Gauss curvature R = a 2 /i 4 . The function h(x) satisfies the equations 
(II. 2p ii) (elliptic type) in the form 

Kikl - „> W x> - cf + O = a lH C + M^coshW* - »)) 
/i x -(/i 2 . - (3/i) 2 (a(x) - c) 2 - 3-M 2 ) = -2-(a/x) 3 -cosh(3/x-(:c - 6)) 

loii/i arbitrary constant C in the first equation; 

(2) h{x) = a-cosh(/i-(x — b)) +c wzi/i some constants fi > 0, a > 0, c, 6 satisfying R = — a 2 /z 4 . 
In case £/ie Killing vector field ^ zs locally a rotation of the hyperbolic plane of 
constant Gauss curvature R = —a 2 fi A and h(x) satisfies the equation (ll.2p ii) ( elliptic 
type) in the form 

M*J - - c) 2 + C) = a,.(C - M*),inh(,,(* - 6)) 

/vOx - (3/i) 2 (a(x) - c) 2 + 3-(/ia) 2 ) = -2-(a/i) 3 -sinh(3/i-(x - 6)) 

m£/i arbitrary constant C in the first equation; 

(3) h{x) = a-sin(//-(x — 6)) + c u>zt/i some constants \i > 0, a > 0, c, 6 satisfying R = —a 2 fj, 4 . 
In this case the Killing vector field zs locally a translation on the hyperbolic plane of 
constant Gauss curvature R = — a 2 /i 4 and a (re) satisfies the equation (ll.2p i) (hyperbolic 
type ) in the form 

h x -(hl + n 2 (h(x) - c) 2 + C) = aa-(C + (a/i) 2 )-cos(,u-(:r - 6)) 

(3.3) 

h x \h\ + (3/i) 2 (/i(x) - c) 2 - 3-(/ia) 2 ) = 2-(a y u) 3 -cos(3 y u-(x - 6)) 

w«£/i arbitrary constant C in the first equation; 

(4) h(x) = a-(x — b) 2 + c some constants a > 0,c,b satisfying R = —4a 2 . In this 
case the Killing vector field is a loxodromy on the hyperbolic plane of constant Gauss 
curvature R = —4a 2 and h{x) satisfies the equation (1 1 . 2 p Hi) (parabolic /nilpotent type) 

h x -{h 2 x - A-a-h(x) + A 2 ) = 2a-(A 2 - 4-a-c)-(x - 6) (3.4) 

with arbitrary constant A 2 ; 

(5) h(x) = a-exp(/ix) + c with some constants fi > 0,a > 0, c and R = 0. In £/ws case the 
Killing vector field -jL is locally a rotation of the Euclidean plane (R = 0) and h(x) 
satisfies the equation (II. 2p i) (hyperbolic type) in the form 

h x -{h 2 x — /i 2 (/i(x) — c) 2 + C) = —afiC-exp(fix) 

h x -(h 2 x - (3/i) 2 (a(x) - c) 2 ) = -8-(a/i) 3 -exp(3/ix) 

with arbitrary constant C in the first equation; 
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(6) h(x) = a-x + c with some constants a > 0, c. In this case R = 0, the Killing vector 
field is locally a translation of the Euclidean plane (R = 0), and h(x) satisfies the 
equation (11 .2p Hi) (parabolic/nilpotent type) 

h x -(h 2 x - A v h(x) + A 2 ) = -a 2 A v x + a-(a 2 - c-A x + A 2 ) (3.6) 

with arbitrary constants Ai,A 2 . 

Proof. As we have shown above, if a metric g admits a Killing vector field v, then in ap- 
propriate coordinates g has the form g = h~ 2 (dx 2 + dy 2 ) and the Killing vector field the 
form v = J^. In this case the Gauss curvature R is given by _R = h xxx -h x — h 2 x . Thus 
we are interesting in possible solutions of the ODE h xxx -h x — h xx = R with constant pa- 
rameter R such that h x ^ 0. By direct calculations we see that every function on the list 
items (l)-(6) satisfies the ODE h xxx -h x — h 2 x = R with an appropriate constant R, and the 
theorem claims that the list is complete. In view of the uniqueness of the solution of an 
ODE with the given initial values we must show that every combination of the initial values 
/ := (R, xq, h(xo), h x (xo), h xx (xo)) is realised by one of the solution on the list. Inverting the 
sign of h{x) and x, if needed, we may assume that h x (xo) > and h xx (xo) > 0. 

Let us consider h xxx (xo) = R+ ^ x ^ x ^ ■ If h xxx (xo) = then the data / are realised by an 
appropriate polynomial of degree 2 (case h xx (x ) ^ 0, list item (4) ) or 1 (case h xx (x ) = 0, 
list item (6) ). 

In the case h xxx (x ) < we set \i := ^-h xxx (x )/h x (x ) and a := yj 'fi- 2 hl(x ) + fi- 4 hl x (x ). 
It is not difficult to see that the ODE h xxx -h x — h xx = R admits the solution h(x) = 
a-s\n(fi-(x — b)) + c (list item (3) ) with appropriate parameters b and c satisfying the initial 
conditions J. 

In the remaining case h xxx (xo) > we set /x := \/ h xxx (xo) / h x (xo) and look for the solution 
of the equation h xxx -h x — h 2 xx = R in one of the forms (1), (2), or (5) with appropriate 
parameters a > 0, 6, c. The form (1) is realised in the case h xx (xo) < fi-h x (xo) in which 
R > 0, the form (2) in the case h xx (xo) > fi-h x (xo) in which R > 0, and form (5) in the case 
h xx (xo) = fi-h x (xo) in which R = 0. The needed parameters a > 0, b, c can be found easily. 

It remains to show every function h(x) given by one of the formulas (l)-(6) satisfies one 
of the ODEs (II. 2p with A ^ and determine possible values of the parameters \x and 
Aq, . . . , A4. Due to the condition A 7^ we may assume that A = 1. The key observation 
is that h(x) is (up to a constant) either a trigonometric (case (3) ), or trig-hyperbolic (cases 
(1) and (2) ), or exponential (case (5) ), or a usual monomial (cases (4) and (6) ) and therefore 
the differential expression h x -(h 2 ±jl 2 -h 2 (x) — Ayh{x) + ^2) will be a polynomial of the same 
type, divisible by the monomial h x , for example, Y2j Bje^ x in the exponential case. 

In the cases (1) and (2) we conclude that the right hand side must be of the form 
Bi-cosh(kfi) +i?2-sinh(A;//) with k = 1, 2 or 3 which gives fx = kii. The case k = 2 is excluded 
by the argument that for k ^ ±1 the expression h x -(h 2 ±jl 2 -h 2 (x)) = h x -(h 2 ±k 2 Li 2 -h 2 (x)) is 
a trig-hyperbolic polynomial of degree 3, i.e., containing a term cosh(3/i) or a term sinh(3/i). 
Using the relations sinh(3x) = 4-sinh 3 (x) + 3-sinh(x) and cosh(3x) = 4-cosh 3 (x) — 3-cosh(x) 
we conclude that the only possible equations are ( 13. ip and ( 13.21) . 

The remaining case (3)-(6) involving trigonometric polynomials, exponential polynomials, 
and usual polynomials instead of trig-hyperbolic ones are treated in the same manner. □ 
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Corollary 3.1. Every equation ( 11. 2ft with Aq 7^ and (A3, Aj) 7^ (0,0) admits only finitely 
many (real) solutions h(x) such that the metric g = h~ 2 (dx 2 + dy 2 ) has constant Gauss 
curvature, except the case of the equation h x (Ao-h 2 + A2) = A4 which always admits a 
solution of the form h{x) = a-x + c with arbitrary c and a satisfying a(A Q -a 2 + A 2 ) = A 4 . 

Every real solution h(x) of fll.2jl is completely determined by its initial values h(xo), h x (xo) 
at a given point xq. Thus for a generic choice of the initial value h(xo) the solution h(x) 
of the equation (11.2j) with this initial value and with any root H x {xq) of the corresponding 
algebraic equation at Xq the metric g = h~ 2 {dx 2 + dy 2 ) has non-constant Gauss curvature. 

Proof. As we have seen, a metric of the form g = h~ 2 (dx 2 + dy 2 ) has constant curvature if 
and only if h(x) is one of the forms (l)-(6). Let us consider possible right hand sides. 

Every expression A 3 - sm ^' x ^ + A^-cos{^-x) can be written in the form A- cos( fj,- (x — b)) with 
unique A and b unique up to a multiple of the period. Similarly, every expression j\ 3 . s ' nh (v-x) _|_ 
A4-cosh(/i-x) can be uniquely written in one of the following forms: 
A-cosh(/i-(x — b)), A-sinh(/i-(x — b)), A-exp(^-x), or A-exp(-^-x). The latter case can be 
reduced to the previous one by inverting the x-axis. Thus the right hand side of the equa- 
tion (ll.2p . ii), determines which type (1), (2), or (5) of the solution h(x) we obtain, and in 
the case (jl.2p . i), the solution must be of the type (3). 

In the case when h(x) is a solution of the type (1), (2), or (3) we proceed as follows: 
Comparing the right hand side of equations ( 11.21) and A3. ID — (13.31) we determine b and possible 
values of /i, there are only finitely many such possibilities. Then multiplying the equation 
by a constant we make A = 1. Next, we compare the l.h.s. and determines the parameters 
c and C . After this the right hand side of (I3.ip - (l3.3p determines the possible values of a. 
Clearly, we have only finitely many possibilities. 

Notice that the type (5) is not generic itself since it occurs only if A 3 = ±fi-A^. Neverthe- 
less, in this case for a given Aq, . . . , A4 we still have only finitely many solutions h(x) giving 
constant curvature. Indeed, we determine possible values of \i considering the right hand 
side of the equation, then from l.h.s. we determines possible values of the parameters c and 
C, and finally again from the right hand side we determine a. 

In the case (4) when h(x) = a-(x — b) 2 + c we must have b = —A4/A3 and a = A\/4, and 
finally c = 2aA g a 2 A:> ■ So for given A , . . . , A4 we could have at most one solution of type (4). 

Finally, if A\ 7^ and h(x) is of type (6), i.e. h(x) = a-x + c, then a must satisfy 
A3 = —a 2 -A\ which gives us at most two possibilities. For every a we have the unique 
possibility for c. □ 

4. Uniqueness of the Principle equation 

The uniqueness of the Principle equation is an interesting phenomenon per se and plays 
an important role in the proof of the main theorem. We shall need the following two results. 

Lemma 4.1. Let h(x) be a complex-valued solution of the equation 

£ ■= h x (h 2 x - 9-h(x) 2 + A 2 ) - A + e 3x - A_e~ 3x = (4.1) 

with complex coefficients A + 7^ 0, A_, A 2 defined for x G [x*, +00). Assume that for x — > +00 
the function h(x) has the asymptotic growth h(x) = a-e x +o(e x ) with a 7^ 0. Then A + = —8a 3 
and there exists a complex-valued real-analytic function f{r) defined for sufficiently small r 
such that f(0) = a and h(x) = e x -f(e~ 2x ). 
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Proof. Write h x (x) = ip(x)-e x , substitute this expression in ( 14. 1 j) . and consider the obtained 
relation as a cubic algebraic equation on a variable ip depending on the parameter x. Then 
for x — > +oo (the coefficients of) the obtained equation converge to ip(if) 2 — 9a 2 ) = A + . 
This implies the asymptotic growth h x (x) = a'-e x + o(e x ) with some a' satisfying the equa- 
tion a' (a' 2 — 9a 2 ) = A + . Integrating it we obtain the asymptotic h(x) = a'-e x + o(e x ). 
Consequently, a' = a and hence a satisfies A + = —8a 3 . 

Now make the substitution x = — |log(r) and h(x) = e x -[a + aoe~ 2x + e~ 4x f (e~ 2x )) = 
r _1 / 2 (a + oqt + t 2 /(r)). Then the equation ( 14. ID transforms into 

8/ 3 r 8 + 3Qff 2 r 7 + 12/ T (3/ 2 + / T a ) r 6 - 12a/ 2 r 5 

-12(3a / 2 + 6// T a + a 2 / r )r 4 + (2A 2 / T - 48aa / T - 72a/ 2 - 36a 2 /)r 3 

+ (3A 2 / + A_ - 72aa f - 8a 3 - \2a 2 f T )r 2 

+ (A 2 — 12aa )(a r — a) = 0. 

This means that we are now looking for solutions f(r) of (14.21) defined for small r > 0. The 
condition on the growth of h(x) and h x means that f(r) = o(r~ 3 / 2 ) and f T (r) = o(r^ 2 ). 
Therefore we must have a = 4p- and the substitution A 2 = 12aa transforms the equation 
( 1421) into 

8/ T 3 r 6 + 36// 2 r 5 + 12/ T (3/ 2 + f T a ) r 4 - 12a/ 2 r 3 

-12(3a / 2 + 6// r a + a 2 / r )r 2 - 12(2aa / T + 6a/ 2 + 3a 2 /)r (4.3) 

+ (A_ - 12a 2 f T - 36aa / - 8a 3 ) = 0. 

For any given A_, a, do, f and sufficiently small r the latter relation can be resolved in f T 
as a real-analytic function f T = F(r, /, a, a , A_) with F(t, /, a, ao, A_) = A ~ 36 "° ^ T ^ 8a o _|_ 
0(r). 

Consequently, we can conclude the following properties: Any solution h(x) of (14.11) satis- 
fying the hypotheses of the lemma is given by the series h(x) = ae x + J^fcLo «fce _ ^ 2fc+1 ^ which 
converges for 16 [xq, +oo). Moreover, the coefficients a, a satisfy the relations 8a 3 = — A + 
and 12aao = A 2 . 

Furthermore, we can also conclude the following two existence results for solutions of ( 14. lft : 

• First, for a, ao satisfying the conditions above, for any given sufficiently large Xq 3> 
and any sufficiently small b G C there exists a unique solution of ( 14.11) with the initial value 
h(x ) = a X0 + a e~ x ° + be~ 3xo . 

• Second, for a, ao satisfying the conditions above and any given a\ there exists a unique 
solution of (14.11) which is defined for x ^> and whose initial terms in the series above have 
coefficients a, a , a±. □ 

Lemma 4.2. Assume that a complex-valued function h(x) satisfies the equation 

h x (A (h 2 x - i?h(x) 2 ) - Ax h(x) + A 2 ) = 0. 

with complex parameters fi, Aq, A\, A 2 . Then R := h xxx h x — h 2 x is constant. 

Proof. We have obviously two possibilities: Either h x vanishes identically, h x = 0, or 

A (h 2 x - /2 2 h(x) 2 ) - A x h(x) + A 2 = 0. (4.4) 

The first case h x = is trivial since then R = h xxx h x — h 2 xx = 0. So we may assume that h x 
is not vanishing identically. 
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Assume additionally that \i ^ 0. In the case Aq = the solution of ( 14.4ft is a constant 
function, and i? vanishes also. In the case A\ = A 2 = the solution of (14. 4p is h(x) = 
C-exp(±fix), and again R vanishes. In the remaining case Aq ^ ^ A\ every solution of 
(14. 4p has the form h(x) = cq + cisinh(±/ix + c 2 ) with arbitrary c 2 and appropriate co,ci. 
This time R must be constant too. 

In the case /1 = the argumentation is changed as follows. If Aq = 0, then h(x) must be 
constant, and then R = 0. If Ai = 0, then must be linear, which is also a contradiction. 
Finally, in the case A ^ ^ Ai every solution h(x) of (14. 4p with \i = is quadratic in x, 
and then R := h xxx -h x — h 2 xx is constant again. □ 

Theorem 4.1. Let h(x) be a complex function defined in a some open set U C C satisfies two 
equations each of the form (11.21) with some complex parameters A , . . . , A4, /1 and respectively 
Bq, . . . , B4, A. Assume that R := h xxx -h x — h 2 xx is not constant. Then \i = ±A and the 
equations are proportional. 

Remark 4.1. For convenience in calculation below we consider only equations of the form 
(II. 2p (ii) or (iii), but not (II. 2p (i). This is an equivalent problem, since the substitution 
fi h-> i/i switches between forms (II. 2p (i) and (ii). 

Proof. One of the techniques to prove the theorem is to write a Taylor series h(x) = 
J2j a j( x ~~ x o) J ') substitute it in both equations, write the expansions, and then compare 
term by term the coefficients. In some places we use another approach, namely, we study 
geometric properties of the solution h(x) using methods of geometric function theory and 
algebraic geometry. It should be noticed however that every relation which will be obtained 
by geometric methods can be also received purely algebraically from the equations obtained 
from Taylor series. 

Denote by E\ and £ M the equations from the hypotheses of the theorem, and A At (x), A^(x) 
the r.h.s.-s of these equations. In the case fi 7^ and resp. A ^ we rewrite them as 
A M = A + e^ x + A_e~^ x and resp. A A = B + e Xx + B_e~ Xx . Notice that by Lemma HH A A (x) 
and A M are non-zero. 

We claim that the function h(x) extends to an analytic multi-sheeted multi-valued) 
function of x G C\5* for some discrete set S C C. Let us consider several cases. The first is 
when Aq 7^ 0. In this case the equation £ M is a polynomial of degree 3 in h x . Let V^x) be its 
discriminant with respect to h x . Then V^{x) = —27A 2 — 4(A fi 2 h(x) 2 — A\h{x) + A 2 ) 3 , this 
is an analytic function in x defined in the domain of definition U of h(x). In the case when 
T>^{x) is not identically zero we can resolve the equation £^ as an analytic multi-sheeted (•£=>■ 
multi- valued) function of x e U\S where S is the set of zeroes of T>^{x). This transforms 
the equation in the explicit form h x = F(x,h(x)) where F is multi-sheeted function with 
ramifications exactly at zero points of the discriminant T>^. This gives us the claim. 

In the case when the discriminant T)^ vanishes identically h x satisfies also the equation 
3h 2 + Aofj^h^x) 2 — Aih(x) + A 2 = which is the derivative of £ M with respect to h x . Again 
we obtain an explicit equation h x = F(x,h(x)) with analytic multi-sheeted right hand side 
F(x, h(x) with singularities in some discrete subset S. Hence this time also h{x) extends to 
an analytic multi-sheeted function of x G C\S for some discrete set 5cC. Finally, in the 
case Aq = the equation £^ can be resolved as h x = jjz^j^ 1 an d we can conclude the 
claim. 
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Assume that A 7^ 7^ B . Then dividing equations by A or resp. B we reduce the 
general situation to the case A = B = 1. 

First, we prove that under the hypotheses of the theorem we must have the relation 
A = ±/x. Let us assume the contrary, i.e., A 7^ ±fi. We shall consider numerous special cases 
and subcases. 

Without loss of generality we may suppose that |A| ^ In particular, A / 0. Ob- 
serve that under complex affme transformation x H- ax + b the equations £ \, £ ll retain their 
structure only changing the parameters /i, A, A^ B^ in particular, A transforms in -. Conse- 
quently, we may assume that A is real positive, A > 0. Then by assumption on \i we conclude 
< A. Further, recall that A A (x) = £> + -exp(Ax) + £L-exp(— Ax) such that at least one 
constant B + ,B_ is non-zero. Inverting the coordinate x, if needed, we can suppose that 

Assume additionally that ji ^ 0. Then A M (x) = v4 + -exp(yux) + A_-exp(— fix). Consider the 
difference £\ — £ ii . It has the form 

h x -(C h(x) 2 + dhix) + C 2 ) = A 5 (x) (4.5) 

where C = A 2 — /i 2 and A^x) = A A (x) — A M (x). Denote this equation by £$■ Integrating it, 
we obtain an algebraic equation 

§C /i(x) 3 + id/^x) 2 + C 2 h(x) = A s (x) + C 3 (4.6) 

which we denote by £$ and in which C3 is some constant and A^x) = J As(x)dx equals 

A$(x) = ^-■exp(Xx) — -^-exp(— Ax) + ^ t -exp( y ux) — ^--exp(— fix). (4.7) 

Making a translation in x, we can suppose that the real axis 1 e 1 does not contain 
singular points of h(x), and that our open set U C C hits the real axis x 6 1. Then there 
exists a unique extension of the function h(x) over the axis x G M. which satisfies the equation 

(USD. 

1 /3 

Since C = A 2 — /i 2 7^ 0, the function h{x) has asymptotic expansion h(x) ~ (Ha) 
exp(Ax/3) for x — )■ +00. From (14. 5 p we conclude that the derivative h x has asymptotic 
expansion h(x) ~ (if^") •|-exp(Ax/3) for x — > +00. 

Further, assume that A 7^ ±3/x. Then h x —n 2 h(x) 2 has asymptotic expansion ~ C-exp(2Ax/3) 
for x — > +00, which gives the asymptotic expansion ~ C-exp(Ax) for the r.h.s. of £ ii . But 
by our condition |3?(/i)| < A the l.h.s. A M (x) has slower growth. The obtained contradiction 
gives the proof in the case A 7^ ±3fi. 

It remains to consider the case A = ±3/x. However, before this case we notice that the 
consideration above are valid also in the case /1 = 7^ A. Indeed, we obtain the same growth 

1 /3 

asymptotic h{x) ~ (|§a) -|-exp(Ax/3) for x — > +00, and the same contradiction in 
growth of right and left hand sides of £ M . 

Now we consider the case A = ±3/i. Recall that we assume that Aq 7^ 7^ Bq. Changing 
the sign of /i, if needed, we obtain A = +3/i. Rescaling x, we can make A = 3 and \i — 1. 
Adding to h(x) a constant we may assume that B\ = 0. At this point we apply Lemma 
14.11 to h(x) and the equation £\. It gives us the presentation of h(x) as a series h(x) = 
ae x + a§e~ x + aie~ Zx + • • ■ which converges for x ^> such that B + = —8a 3 and 12aao = 
B 2 . Differentiating the series we obtain h x = ae x — a§e~ x — 3aie~ 3x + • ■ ■ . This gives us 
hi — h(x) 2 = (h x — h(x))-(h x + h(x)) = —Aaa — &aaie~ 2x + • • • . In the case A\ 7^ the 
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growth of the l.h.s. of is ~ e 2x , which is faster than in the r.h.s. ~ e x . Consequently, 
A l = 0. 

Consider the case i?_ ^ 0. Then repeating the argumentation above, we obtain the 
asymptotic behaviour h(x) ~ a_e _x for x —> — oo. After this we can apply Lemma [4.11 to 
h(x) and the equation S\ in the negative range x G (— oo,Xo], yielding a similar expansion 
h(x) = d-e~ x + a>o-e x + ai~e 3x + • • • , in which £>_ = 8a 3 and 12a„ao,- = B 2 . Under 
translation x H- x + £ with (eC the coefficients a, a_ transform as a \— > ae~^, a_ H- a_e + ^. 
Consequently after an appropriate translation we achieve the equality a = — a_, and hence 
the equality B + = —8a 3 = 8a 3 = The relations B 2 = 12aao = 12a_ao,- gives also the 
equality a _ = — ao- 

The substitution of the series in £ M yields (aA 2 — A + — 4a 2 a )e x + 0(e~ x ') for x — >■ +oo 
and (— a_A 2 — A_ + 4a 2 _ao,-)e _x ' + (^(e 21 ) for x —> — oo. Since vanishes identically, we 
obtain A + = aA 2 — 4a 2 a = —a_A 2 — 4a 2 _a 0> _ = A_. The next term in the expansion of 
for x — > +oo is 

(— 8a 2 ai + 4aQ<2 — a A 2 + 4a 2 a — aA 2 )e~ x . 

Since it must vanish, we obtain ai = — (aA 2 — 4a^a + aoA 2 — 4a 2 ao)/(8a 2 ). Making the same 
computation for x — > — oo we obtain a^- = —a\. 

We conclude that the solution h(x) has sheets which satisfy the relation h(—x) = —h(x) 
Further, recall that h(x) satisfies the algebraic equation £$, see (14. 6p . Since A + = A_ and 
B + = £L, As (given by (14. 7p with A = 3 and /i = 1) is an odd function, A^(— x) = — As(x). 
Consequently, C3 vanishes. 

Let w(z) be the 3-sheeted function of the argument z£C given by the algebraic equation 
|w; 3 + C 2 w = z with C 2 ^ 0. Then h(z) = w(As(x)). Let us observe the following facts 
about the function w(z): The first is that or z small enough three branches of w(z) are given 
by the approximate formulas wq(z) = z/C 2 + 0(z 2 ), w + (z) = +^/—3C 2 /8, and W-(z) = 
— a/— 3C2/8. The second is that the monodromy group of w(z) is the symmetric group 
permuting the sheets of w(z). In particular going along an appropriate path in 2-plane we 
can interchange two given sheets of w(z). 

In turn, the function z = A$(x) is the composition of the polynomial Q(u) := |i? + n 3 + 
2{B + — A + )u with the functions u = sinh(x). The polynomial Q(u) has the same structure 
Cqz 3 + c 2 z as the polynomial |w; 3 + C 2 w . Consequently, there are two possibilities: Either the 
critical values of the polynomials Q{u) and |u> 3 + C 2 w are different, and then monodromy 
of the composition w(Q(u)) in the full symmetric group as in the case of the function w(z), 
or the difference |w 3 + C 2 w — Q(u) splits in the product | n^ifV ~ fyy) with appropriate 
bj G C. In the latter case three possible branches of h(x) are w(u(s\nh(x))) = 6jSinh(a;), 
and in this case R = h xxx h x — h 2 x is constant in contradiction with the assumption of the 
theorem. Consequently, the latter case is impossible, and the monodromy of the composition 
iv (Q(u)) is the full symmetric group. 

Finally, the critical points of the function u = sinh(x) are given by the condition sinh'(x) = 
cosh(x) = 0, and hence the critical values of u = sinh(x) are ±i. The corresponding values 
of z — Q(u(x)) are z± := ±i(|5 + + 2A + ). Here we notice that for each value z± has at 
least two pre-images Q~ l {z±) such that at most one of them is ±i, and neither of these 
pre-images is 0. Further, we notice that the function u = sinh(a;) is surjective. Indeed, it is 
the composition of the rational function u(t) = — t~ l ) and the exponent t = exp(x), the 
map u(t) = \{t — t' 1 ) acts surjectively from C\{0} onto C, the map t = exp(x) surjectively 
from C onto C\{0}. Summing up we conclude that for any of two critical values z* of the 
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polynomial |u> 3 + C 2 w there exists an x* G C which is not a critical point of the function 
z = <5(sinh(x)) but z* = Q(sinh(x*)). This shows that the monodromy of the function h(x) 
is the full symmetric group permuting its 3 sheets(<=> branches). 

In particular, each of 3 sheets of h(x) satisfy both equations and £\. Notice that 
since the l.h.s. A$(x) of H4.6[) vanishes at x — 0, there exists a branch h(x) which vanishes 
at x = 0. Moreover, its behaviour is h(x) = As(x)/C 2 + 0(A 2 (x)). Further, A$(x) = 
2(B + — A + )x + 0(x 3 ) in the case B + 7^ A + , and A$(x) = |£> + x 3 + 0(x 5 ) otherwise. In any 
case, this branch h(x) is regular at x = 0. 

Now consider the equations £ M . By the consideration above, it has the form h x (h 2 — h(x) 2 + 
A 2 ) = 2v4 + cosh(x). This equation has (at most) three local branches of solutions satisfying 
the initial value problem h\ x= Q = each corresponding to a root of the equation 8^=0 
considered as a cubic polynomial on h x \ x= Q. However, we can immediately see that these 
solutions are 2a^sinh(x) where are three roots of the polynomial equation a-(4a 2 + A 2 ) = 
A + ^ The uniqueness of solutions of the initial value problem for ODEs implies the equality 
h(x) = 2a^sinh(x). However, this contradicts to the non-constancy of R = h xxx h x — h xx . So 
finally we have excluded the possibility A = ±3yU with B + 7^ 7^ B-. 

Let us notice that transforming the above solution h(x) = 2asinh(x) by means of afline 
change of the coordinate x {x \- > b(x + xq) with b 7^ 0, xq G C) and the function itself 
(h (->• a'h + c with o'^0,cgC) we can obtain all solutions given in items (1-3) of Theorem 

o 

The next case we consider is when £?_ = 0. Recall that we also have /x = 1 and A = 3. 
We shall consider the asymptotic behaviour of various expressions for x varying some on a 
real line in C given by 3(x) = c and tending to i-c — 00. To simplify notation, we write 
this as x — > —00. Assume that A_ 7^ 0. Then from (14.61) we obtain the asymptotic growth 

h(x) ~ — ic^Y^ e ~ X ^ f° r x ~^ ~ 00 an< ^ f rom (14.51) a similar growth of the derivative h x . 
The substitution of this asymptotic in the l.h.s. of the equation would give the growth 
~ C-e~ x for x — > —00, whereas the l.h.s. £ M decreases. The contradiction shows that we 
must have A^ = 0. Now from (14.61) we conclude that for x — > — 00 the function h(x) is 
given by the converging series h{x) = Yl'jLo a j e '' x w ith some complex coefficients aj. The 
substitution of this series in the equations gives ai(B 2 — 9a>l)e x + 0(e 2x ) = for £\, and 
(A 2 ai — a^ai — A + )e x + 0(e 2x ) = for £ M . In the case a\ = we would have A + = 
and hence A A = 0, which was excluded above. Consequently, B 2 = 9a^. Substituting this 
relation in £\, we consider the further expansion of S\. This gives us aia^e 2x + 0(e 3x ) = 0, 
and hence a = since by the above argument ai 7^ 0. Repeating the substitution we obtain 
-{B + + 8af)e 3x + 0(e Ax ) = from £ x , and {A 2 a x - A + )e x + 0{e 2x ) = for E^. This gives 
us B + = —8a 3 and A + = a\A 2 . Now we see that for each of three roots a± of the equation 
8af + B + = the function a\e x satisfies the ODE £^ and has the correct asymptotic behaviour 
for x — > —00. Consequently, h(x) is one of these three solutions, and hence R = h xxx h x — h 2 x 
must be constant. The obtained contradiction excludes the possibility B_ = 0. 

2 The case when this polynomial has a multiple root is degenerate: In this case the discriminant of 
£^ with respect to h x vanishes identically along the corresponding solution h(x). Moreover, solving the 
corresponding initial value problem in the form of series h(x) = J2j c i x ^ w i tri Co = we obtain the cubic 
equation on Ci-(4cf + A 2 ) = A + which has one double and one simple root. Substituting this double root in 
Ci, all successive equations on c%, C3, . . . can be solved uniquely. This gives us the uniqueness of the problem 
also for this degenerate case. 
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Above we have proven the equality A = [i under hypotheses of the theorem and additional 
assumption A ^ ^ B . Now we consider the case when one of these coefficients vanishes, 
say Bq = 0. Then Bi ^ since otherwise we obtain the equation B h x = A\ whose 
solutions are h(x) = c+e^ + C-e~^ x + c (or h(x) = c 2 x 2 + c\x + Co in the case \i = 0) 
for which R = h xxx h x — h 2 x would be constant. Making transformations h i— >■ ah + c and 
x \-> x/X in the case A^Owe change the equation E\ into h(x)h x = B + e 2x + B_e~ 2x (which 
means that we make A = 2) or respectively h(x)h x = B 3 x + B\. The integration gives 
h(xf = B + e 2x - B_e~ 2x + B 5 or resp. h(x) 2 = B 3 x 2 + 2B 4 x + B 5 with some B 5 E C. As it 
was shown above without loss of generality we can suppose that B + ^ 0. Hence we conclude 
that in the case A = 2 the function h(x) is given by a series h(x) = ae x + Y^jLo a i e ~^ x 
which converges for dt(x) ^ xq ^> and such that a 2 = B + . In the case A = we obtain 
respectively h(x) = ±y/B 3 x 2 + 2B^x + B 5 ~ ±y r B 3 ~x in the case A = 0. In particular, the 

1/2 / 

asymptotic for the derivative is h x ~ _B + e x or respectively h x ~ iyi^. The substitution 
in £^ and comparing of the growth of the left and right hand sides exclude the case A = 
and shows that we could have \i = ±1, fi = ±2, or \i — ±3 in the case A = 2. The case 
/i = ±2 = ±A is our claim, so we must exclude two other possibilities. 

First, we consider the case A = 2 and fi = ±1. As above we distinguish the subcases 
B- ^ and £?_ = and start with the first one B^ ^ 0. Shifting the coordinate x 
appropriately we make B_ = — £>+. Then the coefficients B + , B 5 and a, a , a l5 . . . are related 
as 

R — n 2 R — Onn n — a ' ''~ a ° n — gg(g!-go) n — (a 2 -ag)(« 2 ~5ag) / , a\ 

n + — a iJ 5 — —zaa ai — 2a a 2 — ^ — a 3 — , (4.8J 

and so on. Substitute the series h(x) = ae x + YlJLo Oje - ^- 7 " 1 " 1 ^ and the above relation in 
and write the condition of the vanishing of the resulting expansion. We obtain subsequently 
A + = aA 2 — 4a 2 a , A_ = 4aa^ — 8a 2 ai — A 2 a = 8aa^ — A 2 a^ — 4a 3 , and then the condition 

— — a °^2a 2 12aa °^ = 0. So here we have two possibilities: either A 2 = 12aao or ao = ±a 
(or both). However, after substitution of the first relation A 2 = 12aa in £ M the first non- 
trivial term will be — - e which leads to the relation a = ±a dropped above. But 
then all higher coefficients 01,02,... must vanish and the solution h{x) of the equation 
h(x) 2 = B + e 2x — B_e~ 2x + B 5 will be 2acosh(2x) in the case ao = +a or respectively 
2asinh(2a;) in the case ao = —a. In any case R = h xxx h x — h xx will be constant. The 
obtained contradiction excludes the possibility A = ±2/x, B = 0, and Z?_ 7^ 0. 

Our next subcase is A = 2, \x = 1, and Bo = B- = 0. The procedure here is essentially 
the same as in the previous subcase: Substituting the series h(x) = ae x + YlJLo cije~^ +1 ^ x 
in h(x) 2 = B + e 2x + B§ we obtain the relations 

B + = a 2 B 5 = -2aa ai = - jj. a 2 = ^ a 3 = -jg}, (4.9) 

and so on. Next we substitute the series h(x) = ae x + Y^ C jLo a j e ~ <y2 ^ +lS>x an d the obtained 
relations in and get A + = aA 2 — 4a 2 ao, A_ = 8aa^ — A 2 a , and then the condition 
3a o( A i- 12aa o) _ a s before, setting A 2 = 12aao in £^ we then obtain = which gives 
us the condition a = dropped before. So we must have a = and h(x) = ae x , and hence 
R = h xxx h x — h 2 x will vanish identically. The contradiction excludes also this subcase. 

Next we consider the case case A = 2 and [i = ±3 and start with the subcase one B_ ^ 0. 
As in the case /1 = ±1 above we can additionally assume B- = —B + . Then we obtain the 
same expansion h{x) = ae x + J2T=o a 3 e ~ [2j+1)x with the 

same relations (14. 8p . Substituting 
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them in we obtain subsequently the relations A + = —8a 3 , A\ = 0, A2 = 12aao, A- = 
12a 2 ao— 4oq and then the condition 15 ^° a °) = 0. As above, in both cases ao = ± a all higher 
coefficients a%, a 2 , 03, . . . vanish, the solution h(x) must be either 2acosh(a?) or 2asinh(x), and 
the function R = h xxx h x — h xx will be constant. 

In the subcase A = 2 and u = ±3 and B_ = we obtain respectively first the relations 
( 14. 9p . then subsequently the relations A + = —8a 3 , A\ = 0, A 2 = 12aa , A_ = — 4a 3 ,, and 
then the condition i^ 2 - = 0. The rest follows as in the case A = 2, ji — 1, and B = B_ = 
considered above. 

This finishes the proof of the fact that under the hypotheses of the theorem one has 
the relation u = ±A. Now we show the complete assertion, namely, the uniqueness of the 
equation up to constant factor. As before, we suppose that h(x) satisfies two equations, for 
which we maintain the above notation £ M , S \, A M , A\, Ao, . . . , Bq, . . . , A±, B±. Besides, we 
may assume the equality ji = A. 

Since \i — A, a linear combination of and £\ is again an equation of the same form 
with the same \x. In particular, we can replace £ M or E\ by such a linear combination. 
Consequently, we can assume that Bq = 0, and in the case Aq = Bq = we may also 
suppose that Bi = 0. However, in the latter case we would have B 2 h x = B + e^ x + B_e~ fix 
(resp. B 2 h x = B 3 x + B$ in the case /i — 0) and hence R = h xxx h x — h xx would be constant. 
The contradiction shows that we must have Aq 7^ ^ B±. Normalising, we can make 
A = 1 = Bl 

First, let us consider the case /x = A 7^ 0. Here we apply essentially the same arguments 
as in the above cases A = 2,/i = 1 and A = 2,/i = 3. As we have shown above, making 
appropriate transformations the equation E\ can be brought to the form h(x)h x = B + e 2x + 
B_e~ 2x ) with B + 7^ 0, in particular, we make /i = A = 2. In this way we obtain the algebraic 
equation h(x) 2 = B + e 2x — B_e~ 2x + B 5 and the asymptotic growth h(x) = ae x + 0(e~ x ) 
and h x = ae x + 0(e~ x ) for x — > +00 with a 2 = B + 7^ 0. The substitution gives the growth 
— 3a 3 e 3x + 0(e 2x ) of the l.h.s. of £ M , which contradicts to A M = A + e 2x + A_e~ 2x . 

The argumentation in the case \x = A = is as follows. The equation h(x)h x = B^x + B4 
integrates to h(x) 2 = B%x 2 + 2B^x + B^. An appropriate affine transformation of x and a 
rescaling of h bring this equation into one of the following forms: h(x) 2 = x 2 + 1, h(x) 2 = x 2 , 
h(x) 2 = x, or h(x) 2 = 1. In the cases h(x) 2 = x 2 and h(x) 2 = 1 the expression R = 
h xxx h x — h 2 x vanishes in contradiction to the hypothesis of the theorem. In the remaining 
cases the function h(x) can not satisfy the equation h x (h 2 — A\h{x) + A 2 ) = A%x + A±. 

It remains to consider the case Aq = = Bq ( "Darboux-superintegrable case"). Then both 
A\ and B\ must be non-zero since otherwise R = h xxx h x — h 2 xx would be constant as we have 
shown above. Normalisation of the equations transforms them into h x (h(x) + A 2 ) = A + e^ x + 
A_e~^ x (or = A3X + A4 in the case \x = 0) and respectively h x (h(x) + B 2 ) = B + e Xx + A_e Xx . 
The subsequent integration gives 

^ + A 2 h(x) = ^ x - ^e->" + A 5 , 

^f- + A 2 h(x) = ^fx 2 + A^x + A 5 in the case /i = 0, (4.10) 

+ B 2 h{x) = ^e Xx - ^e- Xx + B 5 . 

In the case A + = A_ = the function h(x) must be constant which contradict the hypotheses 
of the theorem. So one of these coefficient must be non-zero, and changing the sign of /i if 
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needed we can suppose that A + ^ 0. By the same argument B + is non-zero. Observe that 
the equations 14.101 establish an algebraic dependence between the functions e^ x and e Xx in 
the case \x ^ 7^ A, and between the functions x and e Xx in the case \x = 7^ A. This can be 
possible only if \i — ±A. In this situation the difference of the integrated equations 14.101 is 

(A 2 - B 2 )h{x) = *±zl± e ^ _ A^_ e ^ x + (A _ Bg) 

or respectively 

(A 2 - B 2 )h(x) = ^x 2 + (A, - B A )x + (A 5 - B 5 ) 

in the case A = \i = 0. Now it is obvious that the triviality of these relations is the only 
possibility to avoid the contradiction with the condition R = h xxx h x — h xx 7^ const. This 
means the desired proportionality of the equations. 

The theorem is proved. □ 

4.1. Real solutions. Recall that the Principal equations ( 12. 8p . (I2.13P have the following 
meaning: If a surface metric g admits a linear and a non-trivial cubic integral then it has 
the form h~ 2 {dx 2 + dy 2 ) with a function h(x) satisfying one of these two equations with 
some complex parameters fi, Aq, . . . , A4. Of course, we are interested only in solutions for 
which h x is real. In this case h(x) = h\(x) + i-c with some real function h\{x) and a real 
constant c. Substituting we see that h±(x) satisfies the same equation with new parameters 
Aq, . . . , A4. Thus we can consider only real solutions h(x). 

Theorem 4.2. Assume that the equation 1 11. 2ft with some complex parameter fi and complex 
coefficients admits a real-valued solution h(x) such that R = h xxx -h x — h 2 x is non-constant. 
Then \x is real or purely imaginary (or zero) and the equation is complex proportional to 
another equation fl 1 . 2 [) with the same parameter \i and with real coefficients A , . . . , 

Proof. The result follows immediately from Theorem 14.11 applied to the equation (II. 2p and 
its complex conjugate. □ 



5. Number of independent cubic integrals. 
Proof of Kruglikov's "big gap" conjecture. 
Summary of the proof of the main theorem. 

5.1. Number of cubic integrals and Kruglikov's "big gap" conjecture. In |Kr] Krug- 
likov conjectured that the dimension of the space of cubic integrals of a surface metric g of 
non-constant curvature is at most 4. In this section we prove this result for metrics satisfy- 
ing the hypotheses of the main theorem. Our proof applies also for Darboux-superintegrable 
metrics, however, the result in the case is not new. 

Theorem 5.1. Let a function h(x) satisfy one of the equations (II. 2p with complex parameters 
fi, A , . . . , A A . Assume that R := h xxx h x — h 2 x is non-constant. Set H := ^h 2 (p 2 +p 2 ). Then 
the space of complex-valued functions F(x,y;p x ,p y ) that are cubic in momenta (p x ,p y ) and 
satisfy the equation {H,F} = is A- dimensional (as vector space overC). 
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Proof. We distinguish two main cases: fi 7^ and /x = and start with the first one. Set 
L := p y . Then {H, L} = 0. This gives us the following 4 linearly independent solutions of the 
equation {H, F} = 0: L 3 , H-L, and 2-dimensional space of solutions F given by the formulas 
(jl.3p Fl So the theorem claims that there are no more linearly independent solutions. 

We call functions F(x,y;p x ,p y ) satisfying the hypotheses of the theorem (complex) cubic 
integrals (of the Hamiltonian H given by the function h). Denote by the space of complex 
cubic integrals. It was shown by Kruglikov |Krj that the space Th is finite-dimensional@ The 
Jacobi identity implies that the formula L : F i— >■ {L, F} induces a well defined homomor- 
phism C : Th J^h, see §2.11 Consider the decomposition of Th into generalised eigenspaces 
of £ and the corresponding Jordan blocks. Then L 3 and H-L are eigenvectors with eigen- 
value 0. Further, the functions F + and respectively _F_ given by formula (11.31) with C_ = 
and respectively C + = are eigenvectors of C with eigenvalues ±/i. 

It follows immediately from Theorem 14.11 that the space J~h contains no eigenvectors of C 
with eigenvalue A ^ ±/x. Thus in the case /i ^ the assertion of the theorem is equivalent 
to the non-existence of a generalised eigenvector of L with eigenvalue ±/i and the Jordan 
block (Yi)- 

Assume the contrary. Then there would exist cubic integrals -Fo,-^i £ J~h satisfying 
{L, Fx} — ±fiF\ + Fq and {L, Fq} = ±/i-Fo. Inverting the y-axis we can change the sign. So 
we assume that we have +/i in the formulas. Recall that L = p y correspond to the vector 
field J^. Integrating the equations above we obtain F = e m G§ and F\ = ye^ y G + e m G\ 
where G ,Gx are some complex functions of (x,p x ,p y ) cubic in momenta (p x ,p y ) and in- 
dependent of y. Since F is a cubic integral and an eigenvector of C with eigenvalue /i, it 
has the form (11.31) . Since {F ,H} = {H,e m Go\ = 0, the equation {Fx,!!} = now reads 
{e^Gx,H} + {y,H}e^G = 0. Write G x = E'= ^(^"X- Since {y,H} = -p y h 2 x , we 
obtain the equation 

{e^Gx, H} - p y hle^G = (5.1) 
Solving this equation we apply the same procedure as in § 12.21 The bracket {e m Gx, H} 
is given by (12 .4p in which we need to replace a,j(x) by bj(x). Thus the equation (15.11) is 
equivalent to 5 equations which are inhomogeneous versions of 5 equations in (12. 4p with the 
r.h.s.-s given by p y h 2 x e m Go. As in § 12.21 we solve successively the first 3 of them and resolve 
bz(x) from the last one. This gives the following formulas (compare with (12. 5p ): 

bo(x) = B h 3 x 

b 1 (x) = (-(nB + A o yh(*) + %>hl 

b 2 (x) = \-{-{Bx + f)Mx) + (^Bo + 2fxA )-h(x) 2 + 3B h 2 x + B 2 )-h x 
h ^ x ) = i?i 3 - h livBo - A ) - Bx-h(x) + fi 2 -(fiB + A )-h(x) 2 + (fj,B 2 - A 2 ))-h xx , 
Substituting them in the remaining term of (15. ip we obtain the equation 

{3-({iB - A )-h 2 x + /i 2 -(/i-B + A )-h(x) 2 - fi-Bx-h(x) + fi-B 2 - A 2 )-h xxx 
+ e-^o - A )-h x -h 2 xx + (6-/i 2 -(A + /j,-B )-h(x) - 3-fi-Bx)-h x -h xx 
+ Q3-fi 2 -(A + fi-B )-h 3 x 

+ (/i 4 -(/i-5 + 3-A )-h(x) 2 - fi 2 -(fi-Bx + 2-Ax)-h{x) + fi 2 -(fi-B 2 + A 2 ))-h x = 



3 The fact that the parameters pi, Aq, . . . , A4 in Theorem ll.il are real plays no role here. 
4 The proof in [Krj is given for the case H(x,y;p x ,p y ) = (dx 2 + dy 2 )/X(x,y) with real X(x,y). It works in 
our situation without changes. 



TWO-DIMENSIONAL METRICS WITH LINEAR AND CUBIC INTEGRALS 



23 



(5.5) 



which is the counterpart of ( 12. 6ft . As the equation ( 12. 6ft . the above equation can be partially 
integrated in the sense that it can be written in the form (compare with ( 12. 71) ) 

+ (hAh 2 x -B + (^-B + 2^A,)-h{x) 2 - {B l + f)-h(x) + B 2 )) 
+ {£- V 2 ) {h x ihl-A + fi 2 -A -h(x) 2 - A r h(x) + A 2 )) = 0. 
Let us now observe that the other equation {F ,H} = {H, e My G } = is equivalent to 
the equation (h x -(h 2 x -A + fj?-A -h(x) 2 - A v h(x) + A 2 )) = A 3 ^^- + A 4 cos(/ix) with some 
constants ^3,^44, and that the l.h.s. of this equation appears in (15 .4p . Then 

(~cE? ~~ A* 2 ) (^3 Sin (^ x ) + ^4 cos(fix)) = — 2 / u 2 (v4 3 sin(/ix) + A 4 cos(fix)) 
(here we set A' 3 — — ) and so (15 .4p is equivalent to 

h x -(h 2 x -B + (fi 2 -B + 2^A Q )-h(x) 2 - (B 1 + f )-h{x) + B 2 ) = 

(B 3 + y4 4 -x)-sin(/ix) — (A' 3 -x + B 4 )-cos(fj l x) 

For convenience in future let us make the substitution \x H- \fi, Bi 1— )■ i_Bj in the equations 
( 15. 5p and ( 12. 8ft . and rearrange their r.h.s.-s. Then the equations transform into 

h x - (h 2 x -A - n 2 -A -h(x) 2 - A v h{x) + A 2 ) = A + -e» x + A_-e~» x (5.6) 

h x -(hl-B -(fi 2 -B -2^A )-h(x) 2 -(B 1 -f )-h{x)+B 2 ) = 

(B+ + A+-x)-e^ + (B_ - A_x)-e-» x 

Notice that the condition of non-triviality of Fq is equivalent to the non-vanishing of at least 
one parameter A ,Ai,A 2 . Further, by Lemma |4~21 both A + and A_ can not vanish together. 
Replacing fi by — fi, if needed, we can suppose that A + 7^ 0. 

We consider several subcases. The first one is Aq = 0. Then the equation (I5.6P can be 
integrated as 

- 4±h(x) 2 + A 2 h{x) = ^f-e» x - —e'^ + A 5 (5.8) 

If, moreover, A\ = 0, then h(x) = —^-^ x — ~jjx^' e ~^ x + ji^i an d then R = const in this case 
in contradiction with the hypotheses of the theorem. Otherwise we make the substitution 
x 1 — y 2x/jJL. After the substitution ji transforms into 2 and the r.h.s. of (15. 8 p into 4r-e 



2,r 



4p-e 2x + A 5 . So we can conclude that for x — > +00 the function h(x) is given by the 
converging series ae x + Y^=o a i e ■ But then the substitution of this series in 15.71 gives the 
following leading terms for x — > +00: —3a 3 e 3x for the l.h.s., and 2 A±xe 2 ^ for the r.h.s. The 
obtained contradiction shows that the case A Q = is impossible. 

In the case Aq 7^ we make the substitution x t— > 3x/fi which makes /i = 3, and subtract 
( 15 .6p from ( 15. 7p with coefficients B /A Q . This gives us the equation 

h x {3C h{x) 2 + 2C 1 h{x) + C 2 ) = A + xe 3x + C A e 3x + {C 5 x + C 6 )e~ 3x (5.9) 

with some constants Co, Ci, . . . such that Cq 7^ 0. Integrating it, we obtain 

C h(x) 3 + dhix) 2 + C 2 h{x) +C 3 = ^xe 3x + C 4 e 3x + (C 5 x + C & )e~ 3x (5.10) 

with some new constants C 3 , C 5 , C 6 . From this equation we conclude that h(x) is a 3-sheeted 
function on C with a ramification on some discrete subset S C C and that for x —> +00 every 
branch of h(x) has a behaviour h(x) = ax 1 ^ 3 e x {l + A(x)) for some function A(x) admitting 
a converging series Ylij=o dij^~ %x x~'^ 3 with aoo = 0. Moreover, the derivative h x is given by 
the derivative of the expansion above and is a similar series h x = ax 1 ^ 3 e x (l + A(x)) with 
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=o a ij e %x% suc h that a 00 — 0. Substituting these expansions in (EE]) we 
obtain the term — 8a 3 xe 33:; for the l.h.s., which contradicts to the growth A + e 3x of r.h.s. 

This prohibits the possibility A ^ for a solution /z(x) of the pair of equations f l5.6p - fl5.7p . 
and thus excludes cubic integrals Fx, F such that {L, Fi} = F and {L, F } = 0. 

Now we consider the case \i = 0. As above, denote L = p x , set C(F) := {L, F} for any 
function F(x,y;p x ,p y ) and let Th be the space of complex cubic integrals oh H. Then as 
above Th is finite dimensional and £ : Th — > Th is a well-defined homomorphism. It follows 
from Theorem 14.11 that in the case \i = the homomorphism £ : Th — > Th has unique 
eigenvalue \x = and J 7 ^ is a sum of Jordan blocks with eigenvalue \i = 0. 

Notice that L 3 and L-if are eigenvectors of £ with eigenvalue \x = 0. We are going to 
prove that J-^ contains only two linearly independent eigenvectors, a unique Jordan block of 
size 3x3, and no other Jordan blocks. 

Let F Q = F (x, y; p x ,P y ) be given by ( 12. 2p with coefficients a (x), . . . , 03(0;) given by (12.101) 
with Ax = A 3 = 0. Then F is a cubic integral, £(F ) = A V L 3 +A 3 -L-H ^ 0, and C 2 (F ) = 0. 
Assume that we have some other non-zero cubic integral F' such that £ 2 (F') = 0. Then the 
calculation made in § 12.11 and § 12.31 show that F Q must be given by the same (12. 2p with new 
coefficients a' Q (x), . . . , a' 3 (x) given by ( 12.1 0ft with parameters A' , . . . , A' 3 instead of Aq, . . . , A 3 , 
such that h(x) satisfies the Principle equation ( II. 2ft (iii) with A, replaced by A\. 

At this point we obtain two subcases. The first is when both A[, A' 3 vanish, which means 
that JC(F') = 0, i.e., F' is an eigenvector of £ with eigenvalue fj, = 0. In this situation from the 
equation ( 11.21) (iii) we see that either h x is constant or all parameters A' Q , . . . , A' 4 must vanish. 
The first possibility would yield R = h xxx h x — h xx = 0, which contradicts the hypotheses of 
the theorem. Thus all parameters A' , . . . , A' 4 must vanish, and then F' = A[-L 3 + A' 3 -L-H, 
a linear combination of L 3 and L-H. 

Let us underline, that the latter argument demonstrates that the space of eigenvectors of 
£ in Th is 2-dimensional with a basis L 3 , L-H. 

The remaining subcase is when not all parameters A' , . . . , A' A vanish and we obtain a new 
equation of the form (II. 2p (iii). In this situation the uniqueness from Theorem 14.11 ensures 
that A[ = c-Ai with some coefficient c. But in this case F' = c-Fq + A[-L 3 + A' 3 -L-H with the 
same coefficient c and some parameters A[,A 3 . This demonstrates that the space of cubic 
integrals F' satisfying £ 2 (F') = is 3-dimensional with a basis F , L 3 , L-H. In particular, 
we can not have two distinct Jordan blocks. 

Finally, let us show that there does exist a Jordan block of size 3x3, and no Jordan block 
of size 4x4. For this purpose we try to find a cubic integral F satisfying £ A (F) = 0. Since 
the operator £ acts as the derivation in y, the condition £ 4 (F) = means that F is a 
polynomial in y of degree ^ 3. This means that we can write F in the form 

f = EL E U "„(•'•)//'/<:; '/>i (5-ii) 

with some coefficients ay(x). Writing down the equation {F, H} = and considering its 
coefficients at monomials y l p x ~-'p 1 y we obtain 15 ODEs on functions Oij(x) and h(x). We 
solve them subsequently using the conditions h x 7^ 0, h xx 7^ and substituting the results 
in successive equations. Doing so we obtainj^ the following formulas, in which aj J;x denote 
the derivatives of (ii,j(x) and A^ are integration constants: 

(1) a 0>0 {x) = A 00 -h 3 x from 2-h x -(-h x -a 0fi . x + 3-a 0fi (x)-h xx ) = 0; 



'In calculation the authors used Maple® software. 
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(2) a 1)0 (x) = A w -h 3 x from 2-h x -(-h x -a 1)0 . x + 3-h xx -a 1)0 (x)) = 0; 

(3) a 2 , (x) = A 20 -h 3 x from 2-h x -(-h x -a 2>0 - x + 3-h xx -a 2)0 {x)) = 0; 

(4) a 3)0 (a:) = A m -h 3 x from 2-h x -(-h x -a 3>0 - x + 3-h xx -a 3>0 {x)) = 0; 

(5) a 3i i(x) = A 3V h 2 x from 2-h x -(-h x -a 3>1 . x + 2-a 3)1 (x)-h xx ) = 0; 

(6) ai )3 (x) = ^-a 0t2 (x) from 2-/i a .-(-/vai )3 (a;) + h xx -a 0)2 (x)) = 0; 

(7) a 2j3 (x) = ^-a 1>2 (x) from 2-h x -(-2-h x -a 2t3 (x) + /i^-ai^C^)) = 0; 

(8) a 3)3 (x) = |j-a 2!2 (y) from 2-/v(-3-/va 3 , 3 0z) + h xx -a 2>2 (x)) = 0; 

(9) a 3j2 (x) = from 2-h x -h xx -a 3>2 {x) = 0; 

(10) a 0)1 (x) = (-A w -h(x) + A m )-h 2 x from 2-h x -(-h x -A w - h x -a Q ,i ;x + 2-a 0) i(x)-h xx ) = 0; 

(11) a 1>1 {x) = {-2-A 20 -h{x)+A 11 )-h 2 x from 2-/i x .-(-2-^-A 20 - h x -a 1>1]X + 2-ai,i(z)-/i a!3 .) = 0; 

(12) a 2s i(x) = (-3-A 3Q -h(x)+A 21 )-hl from 2-h x -(-3-h x -A 3Q -h x -a 2il . x + 2-a 2tl (x)-h xx ) = 0; 

(13) a 2 , 2 (x) = (-3A 3r /i(x) + §A 20 -/^ + A 22 )-/i x from 2-/i a; -(-3-^-A 3 i + S-Z^/w^o - 
h x -a 2>2 . x + /^-^(x)) = 0; 

(*) A 30 = from Q-h x -h xx -A 30 = at y 3 p 2 x p 2 y , 

(14) a , 2 (x) = (A 20 -/i(x) 2 -Ai r /i(x) + |Aoo-^ + Ao2)-^ from 2-h x -(2-h 3 x -A 20 -h(x) -h 3 x -A n + 
3-h x -h xx -Aoo — h x -ao )2 - x + h xx -ao >2 (x)) = 0; 

(15) ai, 2 (x) = (|A 10 -/i2-2v4 2r /i(x)+A 12 )-/i x . from 2-h x -(2-h 3 x -A 2l -3-h 3 x -h xx -A 10 + h x -a h2 . x - 
h xx -a 1>2 {x)) = 0. 

The latter 4 relations yield certain correction in some formulas above: 

(PJ) a 2 ,i(x) = A 2 r/& 
(H) a 3i0 (x) = 0; 

P) a 1)3 (x) = (A 20 -h(x) 2 - A lv h(x) + \A Q0 -h 2 x + A 02 )-/i xx ; 
dZf) a 2 , 3 (x) = (f A w -h 2 x - A 2V h{x) + \A l2 )-h xx - 
P) a 3!3 (x) = (|A 20 -^ - A 31 -/i(:z) + fA^)-/^. 

Finally, the following formula will be obtained later, we write it here simply for completeness: 

(16) a 0>3 (x) = (-A 10 -h(x) + A m )-h 2 x - \A 23 -x 2 - \A 2i -x + A 03 . 

After calculation of (|1])-( |T5|) and (*) it remains 4 equations. One of them — the coefficient 
at y°p x Py — can be written as 

a , 3 -, x = 2-(A 01 - A 10 -h(x))-h x -h xx - \-A 10 -h 3 x + (2A 2V h(x) - A 12 )-h x (5.12) 

and will be treated later. Three other are the coefficients at y x p x p 3 y , y 2 V x v\ an d y 3 p x Py- After 
normalisation we obtain 

(3A 00 -h 2 x + 2A 20 -h(x) 2 - 2A lv h(x) + 2A 02 )-h xxx + QA 00 -h 2 xx -h x 

+ 6-(2A 20 -h(x) - A u )-h x -h xx + QA 20 -h 3 x + A-{A 22 - 3A 3V h{x))-h x = 

(3A 10 -h 2 x - AA 2V h(x) + 2.1 :,)•// ,,,,, + QA 10 -h 2 xx -h x - \2A 2V h x -h xx = (5.14) 

{3A 20 -h 2 x - 6A 3V h{x) + 2A 22 )-h xxx + 6A 20 -h 2 xx -h x - 18A 3V h x -h xx = (5.15) 

Double integration of the latter two gives 

h x {A 10 -h 2 x - AA 2V h(x) + 2A 12 ) = A 23 x + A 24 , (5.16) 
h x -(A 20 -h 2 x - 6A 3V h{x) + 2A 22 ) = A 33 x + A u , (5.17) 

where A 23 , A 24: , A 33 , A 3i are new integration constants. Substituting the latter relation in 
the equation (I5.12p we can integrate it getting the formula (fl6|) for (1q 3 (x) promised above. 
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The equation (15.131) can not be integrated. However, we can do this even twice with the 
linear combination "( 15.13P " —2- " (15.171) " . and the result is the equation 

h x -(A 00 -h 2 x + 2A 20 -h(x) 2 - 2A lv h(x) + 2A 02 ) = -|A 33 -x 3 - A 3A -x 2 + A l3 -x + A u (5.18) 

with new integration constants y4 13 , A 14 . 

From the above calculation we conclude the following: The metric g = h~ 2 (dx 2 + dy 2 ) 
admits a cubic integral of the form (15. lip if and only if the h(x) satisfies the equations (I5.18p . 
(I5.16p . and (I5.17p . and then the integral can be constructed using the formulas fTTj)- ffl~6]) for 
the coefficients aij(x). 

Observe that the non-zero coefficients a 3j (x) are a 3 ^(x) = A 3 ih 2 and 03,3(0;) = (^A 20 -h 2 — 
A 31 -h(x) + \A 22 )-h xx . This means that the the non-existence of a Jordan block of size 4x4 
is equivalent to the vanishing of all three coefficients A 2 o, A 22 , A31. Suppose for a moment 
that this is not the case. Then the equation (I5.17P is non-trivial. Since we are looking for 
functions h(x) for which R = h xxx h x — h 2 x is non-constant, Lemma T4. 21 says that either A\ 3 
or Au (or both) is non-zero. Consequently, the equation (15.181) is also non-trivial. 

We distinguish several possibilities. The first is when A 20 = 0. Then integrating (I5.17p . 
we obtain 

- 3A 31 -h(x) 2 + 2A 2T h{x) = \A 33 x 2 + A u x + A 35 (5.19) 
with an integration constant A 35 . We can not have A 3 i = since in this case 
would be constant. Also we can not have A 33 = A 3 ^ = by the same reason. In the remaining 
cases we study the behaviour of (branches of) the function h(x) considered as a solution of an 
algebraic equation (15.191) . It follows that for x — > 00 the solution h(x) grows like ~ ax+0(x°) 
(case A 33 ^ 0) or like ~ ax 1 ^ 2 + 0(x°) (case A 33 = 0). From (I5.19P we obtain the growth 
of h x : ~ a + 0(x _1 ) in the case A 33 7^ and ~ ax~ x l 2 + 0(x _1 ) in case A 33 = 0. In any of 
these two cases we see that the growth of the l.h.s. of (15.181) is slower than that of the r.h.s. 
The contradiction shows that we must have A 2 q ^ 0. 

Now subtract the equation (15.171) from (15 . 1 8[) with coefficient We obtain the equation 

h x -(2A 20 -h(x) 2 + A' u -h(x) + A' 02 ) = ~lA 33 -x 3 - A u -x 2 + A[ 3 -x + A' u (5.20) 

with new constants A' u , A' Q2 , A' 13 , A' lA and the same constants A 20 7^ 0, A 33 , A 34 . Integration 
gives 

^■h{x) 3 + ^-h{x) 2 + A' 02 -h{x) = -^f -x 4 - ^f-x 3 + ^f-x 2 + A' u x + A 15 (5.21) 

with an integration constant A^. Our next possibility is A 33 7^ 0. In this case the same 
argument as above gives the asymptotic h(x) ~ ax 4 / 3 + O(x) and h x ~ yx 1 / 3 + O(x ) for 
x — > 00. Moreover, the coefficient a satisfies the relation A 33 = — 8A 20 -a 3 . Substituting these 
asymptotic in (I5.17P gives the asymptotic —8A 3 i-a 2 -x 5 ^ 3 + 0(x 4//3 ) for the l.h.s. provided 
A 3 i 0. This is a contradictions. Consequently, we must have A 3 i = 0. But then we obtain 
the asymptotic ^A 3 i-a 3 -x + 0(x 2 ^ 3 ) for the l.h.s., whereas on the r.h.s. we have the leading 
term A 33 x = — 8A 20 -a 3 x. As the result we conclude that the case A 33 7^ is impossible and 
we must have A 33 = 0. 

Notice, that in this remaining case A 33 = we must have A 34: since otherwise R = 
h xxx h x — h 2 xx would be constant by Lemma I4.2L Here as above from (15. 2 1 1) we can conclude 
the asymptotic h(x) = ax + O(x ) and h x = a + 0(a;~ 1 ) for x — > 00. The substitution of these 
asymptotic in (15.171) gives linear growth — 6A 3 i-a-x + O(x ) for the l.h.s. provided A 3 i 7^ 0, 
whereas the r.h.s. is constantly A 3i . So we must have A 31 = 0. But the only solutions of 
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the equation ( 15. 17ft with A31 = A 33 = are h x = const which contradict the hypothesis 
R = h xxx h x — h xx 7^ const. The latter contradiction demonstrates the non-existence of 
Jordan blocks of size 4x4. 

Finally, we describe cubic integrals F which give Jordan blocks of size 3x3. The latter 
condition means that F is given by (15. lip with vanishing coefficients a^j{x), j = 0, . . . , 3. 
From formulas (TTjl- (fl6il we see that this condition is equivalent to vanishing of parameters 
^20 , ^.31, ^22- In this case we see from ( 15. 17ft that the coefficients A\ 3 , An must also vanish. 
Thus the equation (I5.17P becomes trivial, and the equation (15.181) simplifies to 

h x -(A 0Q -h 2 x - 2A ir h(x) + 2A 02 ) = A 13 -x + A u . (5.22) 

The uniqueness of the equation proved in Theorem 14. 1 1 ensures that the equation (15.161) and 
(I5.22p must be proportional to each other and to the equation (ll.2p (iii). Let us denote by 
Ci,C 2 the corresponding proportionality coefficients. This leads to the following relations: 

Aoo = Ci-Aq, An = h-Ci-Ai, Aq2 = y-CyA 2 , A13 = Ci-Az, An = Ci-A^, 

(5.23) 

A10 = C2'A , A 2 l = |-C2-Al, A% 2 = 2-C2-A 2 , A23 = C2'A 3 , A 2 4 = C2-A/1. 

Additionally set A 03 = Clz and A Q i = \Clh- Substituting these relations and formulas 
(TTT)-(lT6l) we obtain the following formula for a cubic integral F in the case /i = 0: 

F = C L3 -p 3 y + C LH -\hl{plp y + pi) (5.24) 

+ C v (Ao-hl-p 3 x + \-A v hl-plp y + \-(3A Q -h 2 x - A v h(x) + A 2 )-h x -p 2 y -p x 

V N (5-25) 

y r (3A -h 2 x -A v h(x) + A 2 )-h xx -p 3 y j 

+ C r (A Q -hl-ypl + (£-A v h 2 x - A Q -h 2 x -h{x)\p 2 xVy 

+ H3A Q -h 2 x - A r h(x) + A 2 )-h x -yp x p 2 y (5.26) 
+ (li3Ao-h 2 x - A r h(x) + A 2 )-h xx -y 2 - A a -h 2 x -h{x) - \-A A -x - \-A 3 -x 2 )-pl 



Thus F is a linear combination of 4 independent cubic integrals as stated in Theorem 11.11 
two of which are L 3 and L-H, and other two are linear and quadratic in y. 

Notice that the formulas ( I5.25P and (I5.26P for cubic integrals differ from those in (II. 5p and 
(II. 6p . To obtain the last ones we need to subtract from (15.25P and ( 15.26P the derivative of 
the equation ( II. 2p (iii) with coefficients \-p\ and respectively \-p\- □ 

5.2. Summing up: Theorem [TTTJ is proved. As we explained in § §2.1||2.2|I2.3| any metric 
g admitting linear integral L and cubic integral F such that L, H and F are functionally 
independent have in an appropriate coordinate system the form h~ 2 (dx 2 + dy 2 ), where the 
function h(x) satisfies (JZl) or (12~T3|) . 

In Theorem 14.21 we proved that fi in (I2.8P is real, or pure imaginary, and the parame- 
ters A ,...,A4 in ( 12. 8p and in (12 . 131) become real after multiplication with an appropriate 
constant. Thus, the equations ( 12.81 12.131 are essentially the equations ( 11.21) . 

It follows from § §2.1| 12. 2| 12.31 that the metric g = h~ 2 (dx 2 + dy 2 ) admits at least one 
cubic integral i*\ which has the form from Theorem ll.l| and two functionally independent 
integrals F±, F 2 in the case /i 7^ 0. In Section |5] we construct an additional independent cubic 
integral F 2 in the case fi = 0, and explained why there are no other integrals (unless the 
metric has constant curvature). 
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Theorem 11.11 is proved. 

6. Global solutions 

In this section we show that if the function h(x) satisfies the equation (jl.2p (ii) and h'(xo) > 
at some point xq whereas the real parameters fi > 0, Aq, . . . , A4 satisfy inequalities A > 0, 
\i-A± > \A 3 1 then the metric g = -^(dx 2 + dy 2 ) smoothly extends to the sphere S 2 together 
with the linear integral L = p y and the cubic integral F given by (11.41) . More precisely, we 
show that if (r,ip) are polar coordinats on M 2 related to (x, y) by r = e x ^,tp = y/fJL, then 
g, L, and F are well-defined on the punctured plane M 2 \{0} and extend smoothly to the 
origin and to the infinity point 00 of the Riemann compactification S 2 = C = M. 2 U {00}, 
such that the extended tensor g is still a (non-degenerate) Riemannian metric on the whole 
sphere S 2 . 

The obtained in this way family of examples of superintegrable metrics on the sphere is 
new. Indeed, by |Ki] Darboux-superintegrable metrics can not live on a closed manifold, so 
the only known superintegrable metrics on the 2-sphere are the standard metrics of constant 
curvature. In view of Corollary 13.11 for most values of the parameters satisfying the above 
conditions, the metrics are not metrics of constant curvature. 

Let us describe the conditions on parameters fi > 0, A , . . . , A4 which distinguish our 
global solutions. Since we want that the linear integral L to be also globally defined on S 2 , 
the Killing vector field must be as in the standard rotation (see for example |Bo-Ma-Fo] ). and 
we must have the elliptic case , i.e., h must satisfy ( II. 2p . (ii). Since Darboux-superintegrable 
case is impossible on closed surfaces due to |Kij . A is non-zero. Then dividing the equation 
by Aq we obtain Aq = 1. Applying the action h(x) 1— > h(x) + c we can make A\ = 0. Since 
yU ^ 0, making an appropriate rescaling in x we can make \x = 1. Further, we rewrite the 
free term 74 3 sinh(x) + t4 4 cosIi(x) in the form A + e x + A_e~ x and we impose the positivity 
condition A + , A^ > 0. This means that the free term is positive for all 16R, An appropriate 
translations in x direction transforms the term A + e x + A_.e~ x into A e -(e x + e~ x ). It is easy 
to see that in terms of the original parameters our conditions are 

/I >0, A >0, M-n- \M > 0. (6.1) 

Lemma 6.1. For any A e > 0, A 2 and h(xo) there exists a unique real-analytic local solution 
h(x) of the ODE 

S ■= h x (h 2 x - h(x) 2 + A 2 ) - A e -(e x + e~ x ) = (6.2) 

with the initial value h(xo) such that h x (xo) is the unique positive root of the characteristic 
polynomial x(A) := A(A 2 + A 2 - h(x ) 2 ) - A e -(e x ° + e~ x °). 

Proof. Considering the graph of the polynomial A(A 2 + a) for different a (mainly, for the cases 
a ^ and a < 0) we see that there exists a unique positive solution of the equation A(A 2 +a) = 
A with A > which depends real-analytically on a and A > 0. The standard theory of ODEs 
implies now the desired local existence and uniqueness of solutions of (16. 2p . □ 

Remark. Notice that depending on coefficients A > and a one could have two distinct, 
one double, or none negative roots of the equation A (A 2 + a) = A. Vice versa, in the case 
A < there is one negative root and there could be two distinct, one double, or none positive 
roots. This explains our "Ansatz": we need that the right hand side 74 3 sinh(x) + y4 4 cosh(x) = 
A + e x + A_e~ x remains positive. 
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Let h(x) be a local solution constructed in the previous lemma. Make the substitution 
x = log(t). Then the equation ( 16. 2\\ transforms into 

£ := t-ht ({t-h t ) 2 - h(t) 2 + A 2 ) - A e -(t + t~ l ) = (6.3) 

defined for all t > 0. 

Proposition 6.1. Any local solution h(t) of ( 16.31) with h t positive extends to the interval 
t G (0, +oo). Moreover, the functions t-h(t), t 2 h t , and t 3 h a extend real- analytically to some 
interval t G (— e, e), e > 0. In particular, the function h(t) has a simple pole at t = 0. 

Notice that in terms of the variable x the first assertion of the proposition states that 
every local solution h{x) as in Lemma [6.11 is global, i.e., is defined for all 

Proof. Considering the local behaviour of h(t) at t — 0, we show a stronger property, 
namely, that h(t) = f(t 2 )/t for some real-analytic function /(r), defined in some inter- 
val r G [—e, +oo]. This means that we make the substitution t = \/t, or equivalently, t 2 = r 
or x = |logr. However, we consider all three functions h(x), h(t), and h(r) as a single object 
given in different coordinates. 

For A positive and h, A 2 real, denote by r\ = rj(h, A 2 ; A) the unique positive root of the 
polynomial X(X 2 +A 2 — h 2 ) — A. Then rj(h, A 2 ; A) is monotone in every its argument A > 0, A 2 
and h ^ 0: rj(h, A 2 ; A) will increase if we increase A > and the absolute value \h\ and 
decrease A 2 . Consequently, for any given A 2 and A e > any solution h(x) of (16.21) with 
positive h x satisfies h x ^ rf := 77(0, A 2 , ; 2A e ). In particular, h(x) is monotone. 

Now let us observe the following two facts. First, the translation in x (which means the 
multiplication of t by a constant) transforms our equation ( 16. 2h into 

hjhl ~ h{x) 2 + A 2 ) - A + -e x - A_-e~ x = (6.4) 

or respectively (16. 3p into 

t-h t ((t-h t ) 2 - h(t) 2 + A 2 ) - A+-t - A^-r 1 = (6.5) 

with positive parameters A + ,A_. Second, for any values of parameters our equation has 
a very simple function as the solution. Namely, the function h(t) := C/ l (— A^/t + A + -t) 
satisfies the equation (16. 5p if and only if the constant Ch is the unique positive root of the 
polynomial Ch{^A + A_-Cf l + A 2 ) = 1. We use these special solutions and their translations 
in x reparametrisations t 1— > c-t) to estimate the behaviour of a general solution h(t). 

Let h(t) be any solution of (16.31) defined in a neighbourhood of the initial value to- Denote 
h := h\ t=to . Let h{t) be the unique solution of (I6.5P of the form h{t) := Ch(—A_/t + A + -t) 
with Ch > 0. In the case when h(to) = ho we must have h(t) = h(t) everywhere and thus 
h{t) is defined globally on (0, +00). 

Assume that h{to) > ho. Since both functions are solution of the same 1st order ODE, 
h(t) > h(t) for every t from the maximal existence interval (t_, t + ) C (0, +00) of the solution 
h(t). Further, since h(t) is monotone increasing and bounded from above by the function 
h(t) defined on the whole ray (0, +00), we conclude that the solution h{t) does not explode 
and exists on the whole interval (t ,-\-oo). In particular, the existence interval for h(t) is 
(*_, +00). Moreover, h(t) ^ C h A + t for all t G (t_, +00). 

Our next step is construction of a similar lower bounding function /i_ (t) which also has the 
form h_(t) = Ch-(—A'_/t + A + -t) and satisfies the equation (16.51) with some new parameters 
A'_ and A 2 . This means that constructing from h(t) we change only one parameter, 
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namely A_. Increasing it the value h-(to) will decrease. Thus we can find the new value 
A'_ > A- from the condition h^(to) = ho. Then we find the new value A' 2 from the relation 
Ch{4:A + A'_-C^ + A' 2 ) = 1. Notice that since A'_ was increased, A' 2 is lower than A2, i.e., 
A' 2 < A 2 . Assume additionally that ho ^ 0. Then h(t) < for every t G (to, to). Now 
using the monotonicity on the function rj(h, A 2 ] A) and comparing the equations for h(t) 
and h-(t) and the initial values h-(to) = h(t ) and 4/i_(to) > ^.^(to), we conclude the 
following: for every t in the whole existence interval (i_,io) left from the initial point to we 
have h-{t) < h{t), \h_(t)\ 2 > \h(t)\ 2 , and f t h.(t) > f f h{t). 

The mononoticity argument above can be applied in the case when ho > 0. In this case 
we need an additional step. Namely, since we have the uniform estimate h'(t) ^ 77* > 0, 
our solution h(t) vanishes at the unique t\ lying in the existence interval (£_,£o) left from to- 
Then we apply the same argument above at the point t± instead of to. 

Finally, recall that above we have proceeded under assumption h(to) > ho- The remaining 
case h(to) < h is treated quite similarly, and we only indicates the changes. In this new 
situation the function h(t) := C^(— A_/t + A + -t) will be a lower bound for our solution h(t), 
i.e., h(t) > h(t) everywhere in the existence interval. This fact together with monotonicity 
of h(t) will imply the extensibility of h(t) left until the value = 0. To construct the upper 
bounding function h + (t) we increase the parameter A + and also decrease A 2 . Moreover, 
in the case ho < we need an additional step, in which we go right to the unique point 
t\ G (to,t + ) such that h(t\) = 0, and apply the argument at this point t\. 

For the second assertion of the proposition in the case h(t ) < h we need one bounding 
function, namely, the function h + (t) of the familiar form h + (t) := C/ l (— A'_/t + A + -i) such 
that h + (t) > h(t) for t < t Q (h + (t) does this for t > t Q or respectively for t > t\ > t ). As 
before in the case h > we go left to some point t 2 < t such that hfa) < 0. Notice that 
h{t,2) is still less than h(t2), and hence hfa) < 0. Now, if we start to decrease v4_, the value 
h(t2) will increase until it arrives to 0. Consequently, for some A' Q G (0, A ) the function 
h + (t) := Ch{— A'_/t + A + -t) has value h+fa) = hfa). The same monotonicity argument as 
above gives us h{t) < h + (t) for t G (0,t 2 ). 

This gives the desired global existence of the solution h: the maximal existence interval is 
(0, +00) for the coordinate t which means that h(x) exists for all 

As a consequence of the above argument, we obtain the estimate — j < h{t) < — | for 
< t 1 with some positive constants c, C for any solution of (16. 3p . 

Now let us make the substitution t = y/r and h{t) = f(t 2 )/t = f(r)/y/r. Then the 
equation ( 16. 31) transforms into 

8-/ r 3 -r 2 + (-12-/ r 2 -/(r) -A e + 2-f T -A 2 )-r - f(r)-A 2 -A e + A-f(r) 2 -f T = (6.6) 

For r small enough and / < we can resolve this equation with respect to the derivative f T , 
and then (16.61) transforms into the form 

f T = $(Tj(T);A 2 ,A e ) (6.7) 

for some real-analytic function $(r, /; A 2 , A e ) of arguments r G [—£,£], / < 0, A 2 ,A e . 
In particular, $| T= o = A£ \fi ^ ■ It follows, that for every A 2 ,A e , every negative fo, and 
every To small enough there exists the unique solution of the equation (I6.6P with the initial 
value fijo) = fo- More precisely, for given intervals A 2 ,A e G [— C, C], f G [— C, — c] with 
< c < C there exists an e = e(c,C) > such that the equation (16. 6p has the unique 
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solution with the initial value /(to) = fo at r G [—£,£] such that / T (ro) = $(to, /o; ^2, ^4 e )- 
More over, this solution is well-defined and real-analytic on the whole interval r G [— e, e]. 

Next, we observe that the above estimate — j < h{t) < — | for < £ <C 1 is equivalent to 
the estimate — C < /(r) < — c for < t <C 1. The proposition follows. □ 

Now we state the result about the extensibility of g and F to a metric and a cubic integral 
defined globally on S 2 . 

Theorem 6.1. For any A e > 0, A2 and ho, let h(t) be the unique solution of the equation 
( 16. 3 p with the initial value h\ t= i = ho and with ht\t=i > 0. Then the metric 



dt 2 + t 2 -dp 2 

9 = ¥h 2 (6 - 8) 

defined on the 'plane M 2 with the polar coordinates (t, tp) extends to a real analytic metric on 
the sphere >S 2 = C = lR 2 U{oo} with the Killing vector v = which admits a cubic integral 
F also well-defined and real-analytic globally on S 2 . 

Moreover, the metric g has constant curvature if and only if ho = 0. 

Proof. By Proposition 16. II the function h(t) is well-defined for all t G (0, +00) and for t small 
enough h(t) = with some real-analytic function /(r) such that /(0) < 0. Consequently, 
in a neighbourhood of the origin the function t 2 h t is real-analytic and non- vanishing. It 
follows that the formula g = dt defines a non-degenerate real-analytic Riemannian 

metric in a neighbourhood of the origin in M. 2 with the polar coordinates (t, tp). Substitution 
t = e x , (p = y transforms this metric into the familiar form g = dx ^ dy and the equation 

(16. 3 p into (16. 21) . By Proposition 16.11 the metric g = dx X dy is well-defined for all x G M or 

222 

equivalently for all t G (0, +00). This means that the metric g = dt ^ h i dip is well-defined on 
the whole plane R 2 . 

To show the extensibility to the infinity point 00 we apply the inversion of the sphere 
5* 2 = IR 2 U {00} with respect to the unit circle given by the condition t — 1. Recall that 
the inversion map interchanges the origin and the infinity point 00 and that in the polar 
coordinates it is given by (t,<p) H- (t' 1 ,^). Changing to the coordinates x = \og(t),y = ip 
we obtain the formula (x,y) 1— > (—x,y). So we conclude immediately that the extensibility 
of the metric g to the infinity 00 is equivalent to the extensibility to the origin 0, and thus 
this is the case. 

Let £ := t-cos(p) and 77 := t-sin(^) be the Cartesian coordinates corresponding to the 
polar coordinates (t,<p) and p^Pr, the corresponding momenta, i.e., dual coordinates on 
T*S 2 . Then the vector field d y = d v is given by ^d Tj — rjd^. This means that the linear 
integral L = p y is given by L = C,p v — rjp^ and hence L extends smoothly to the origin. By 
the symmetry argument L extends also to the infinity point 00. 

It remains to show that the cubic integral F given by (II ,4p also extends to the origin and 
to the infinity point 00. Our argumentation is as follows. First, we substitute in the formulas 
(II. 4p our values of parameters A = /x = 1 and A\ = 0. Next, without loss of generality we 
may set = for "phase parameter" in (11.41) . Then each function ao(x), . . . , a^(x) in (11.41) 
becomes a sum of linear and cubic monomials in functions h(x), h x , h xx . Now we observe that 
each function h(x), h x , h xx , considered as a function of the variable t, has the form fi(t 2 )/t 
for certain real- analytic function fi(r), namely, /o( r ) = /(t) as in the proof of Proposition 
16. 1\ /i(t) = 2r/'(r) — /(r) and f^ij) = 4r 2 /"(r) + f(r). Further, in the same way as it 
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was done for p y = — rjd^, we obtain the formula p x = + rjp^. Substituting all these 
relations and also the relations cos(y) = |, sin(y) = ^ in fll.4p we see that in coordinates 
(£, 77) the integral F has the form 

where r = t 2 = C, 2 + rf and ipij{r) are some real-analytic functions. It follows that to 
establish the real-analyticity of F we need to know only the linear parts of the functions 
ifjijij). Considering the expressions of these linear parts of ipijij) in terms of the the functions 
h(x), h x , h xx we see that only two lower monomials of each function h(x), h x , h xx are involved. 
We write h{t) = —c -t~ l +c 1 -t+0(t 3 ). Differentiation yields h x = t-h t = +c -t~ 1 +ci-t+0(t 3 ) 
and h xx = t-4f(t-h t ) = — Co-t" 1 + cyt + 0(t 3 ). Substituting these expressions in the equation 
(16. 3p and considering the coefficients by t , we obtain algebraic equations 

4-cg-ci + c -A 2 - A e = 0, A-cj-co + c r A 2 - A e = 0. 

So we conclude the equality C\ = c and the formula A e = c -(4-Cq + A 2 ). Now, substituting 
all these formulas in (11.41) we obtain 

F = cl- Pv (p 2 +p 2 ) + (3-cU 2 + C ° iA T Cl ] - V 2 )-p! - coi2-c 2 + A 2 ).H.p 2 . Pv 

+ r (1 °-f M) -e + eg v)-^-Pe + 2-cg.17.e-pj + o(e + v A ) 

in which 0(£ 4 + t] 4 ) means term of higher degree in £, 77. 

This shows that the cubic integral F extends real-analytically to the origin as desired. 
The extensibility of F to the infinity point 00 can be obtained from the extensibility to the 
origin by means of the inversion. 

The theorem follows. □ 

7. Conclusion 

We found all two-dimensional Riemannian metrics whose geodesic flows admit one integral 
linear in momenta (L) and one integral cubic in momenta (F) such that L, F and the 
Hamiltonian H of the geodesic flow are functionally independent. Within these metrics, 
we point out the metrics that are already known, and proved that most of our metrics are 
new. We have also showed that, in the case when the parameters satisfy certain inequalities, 
the metric and the integrals L and F extend real-analytically to the sphere S 2 , giving new 
unexpected examples of integrable metrics on the sphere. 

The results and the methods of our paper suggest the following directions of further 
investigations. 

Problem 1. Generalise our result for integrals of higher degree. 

In other words, we suggest to construct all two-dimensional metrics whose geodesic flows 
admit one integral linear in momenta (L) and one integral polynomial in momenta of degree 
4, (5, 6, etc.) in momenta (F) such that L, F and H are functionally independent. 

The main trick that allowed us to solve the case (linear integral L + cubic integral F) 
survives in this setup: the Poisson bracket {L, F} is again an integral of the same degree 
as F . Arguing as in §2.11 one can reduce the problem to analyse of certain system of ODE. 
Though it is not clear in advance whether one can reduce this system of ODE to one equation 
(as we did in the case (linear integral L + cubic integral F)), the approach should at least 
allow to construct new examples of superintegrable metrics. 
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Problem 2. Generalize our results for pseudo-Riemannian metrics. 

We expect that this is possible to do the local description using the same idea. We 
do not expect that one can find the analog of our global examples on closed surfaces in the 
pseudo-Riemannian case. Generally, it could be complicated to generalize global Riemannian 
construction to the pseudo-Riemannian setting. In certain cases though the existence of 
additional structure such as additional integrals (as for example in |Ma3j ) allows to keep 
control over the situation. 

Problem 3. Quantise the cubic integral. 

Take a metric g from Theorem 11.11 and consider its Laplacian A (since our metric is 

^{dx 2 + dy 2 ), A = h 2 x \Jf^ + ; one can view it as a mapping A : C°°(M 2 ) -> C°°(M 2 ) 

though one also can consider Laplacian as a linear operator on bigger function spaces). 
Does there exist a differential operator F of degree 3, 

F = a {x,y)^ + a 1 {x,y)^^ + a 2 {x,y)£^ + a 3 {x,y)-§^ 

+ b (x,y)^ + b 1 (x,y)£-j^ + b 2 (x 1 y)^ + c Q (x,y)£ + c 1 (x,y)£ + d(x,y) 

such that F commute with A, i.e., such that for every smooth function / : M 2 — > R 

[A,F](f) := A(F(/)) - F(A(f)) = 0, 

and such that its symbol 

«o(z, y)pl + ai(x, y)v\v y + av{x, y)p x p\ + a 3 (z, y)p 3 y 

coincides with the integral F from Theorem 11.11 ? 

Recall that for all previously known superintegrable systems such quantisation of the 
integral was possible, and was extremely useful for the describing of eigenfunctions of A. 
Note that quantum superintegrability can be most effectively used (if it exists) in the case 
of metrics from Theorem 16.11 since in this case the Laplacian is a selfadjoint operator. 

Problem 4. Find physical or mechanical systems realizing the Hamiltonian systems corre- 
sponding to (at least some) metrics constructed in Theorem \l.l\ 

This problem is an interesting challenge for both mathematicians and physicists, especially 
in the case of global systems given by Theorem 16.11 Let us note that many classical examples 
of global integrable systems have arisen as the mathematical models for concrete naturally 
defined dynamical systems in physics and mechanics, and that many superintegrable metrics 
have physical realisation or can be applied to solving of physical problems. 

Problem 5. Describe the metrics from Theorem 1 6. 1\ in the terms of [Bel Chapter 4]. 

As we already mentioned above, all geodesies of the metrics from Theorem 16. II are closed. 
By construction, the metrics are the metrics of revolution. Then, these metrics are a subclass 
of the so-called Tannery metrics from [Bel Chapter 4]. 

Problem 6. Find isometric imbeddings of metrics from Theorem 1 6. 1\ in (M. 3 ,g stan dard)- 

Such isometric imbeddings are possible at least for certain metrics from 16.1 1 since they have 
positive curvature as small perturbations of the standard metric. This problem is related to 
Problem [U and is a geometric analog of it. It also is related to Problem [5] in view of |Be[ 
Chapter 4(C)]. 
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